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ABSTRACT 


An investigation into fiber fracture and debonding in metal matrix 
composites is conducted using the finite element method. The superelement 
finite element technique was used to model a metal matrix composite under 
various loading condition and with varying degrees of fiber debonding. The use 
of superelements saved many man hours by allowing for alteration of only the 
primary superelement to manipulate partial bonding for the entire model. The 
composite’s material properties were calculated and the effects of fiber debonding 
on these properties were noted. The internal stress state of the composite while 
under various loads was also investigated. Special interest was devoted to the 
change in stress state as a result of increasing fiber debonding. 
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I. 39% Debonding, cr e Loading 

yy 

A.66 a Normalized Microstresses, O-C Direction, 

li ’ 

2.78% Debonding, a e Loading 
yy 

A.67 a Normalized Microstresses, O-C Direction, 
n 

4.17% Debonding, cr e Loading 
yy 

A. 68 <r Normalized Microstresses, 0— C Direction, 

li ’ * 

5.56% Debonding, <r e Loading 
yy 

A. 69 a Normalized Microstresses, 0— C Direction, 

li * ’ 

6.94% Debonding, <r e Loading 
yy 

A. 70 a Normalized Microstresses, 0— C Direction, 

li 

8.33% Debonding, <r e Loading 
yy 

A.71 a Normalized Microstresses, 0— C Direction, 

n ’ 

9.72% Debonding, a e Loading 
yy 

A.72 a Normalized Microstresses, 0— C Direction, 

li ’ 5 

II. 11% Debonding, a e Loading 

yy 

A. 73 a Normalized Microstresses, 0-B Direction, 

0.0% Debonding, a e Loading 
yy 

A. 74 a Normalized Microstresses, 0— B Direction, 

1.39% Debonding, a* Loading 
yy 
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255 
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XXIV 



A.75 g Normalized Microstresses. 0— B Direction, 

33 1 

2.78% Debonding, <r e Loading 261 

yy 

A. 76 g Normalized Microstresses. 0— B Direction, 

33 

4.17% Debonding, <r e Loading 262 

yy 

A. 77 g Normalized Microstresses, 0-B Direction, 

33 ’ 

5.56% Debonding, g e Loading 263 

yy 

A. 78 g Normalized Microstresses, 0— B Direction, 

33 1 

6.94% Debonding, cr e Loading 264 

yy 

A.79 g Normalized Microstresses, 0— B Direction, 

33 ’ 

8.33% Debonding, G e Loading 265 

yy 

A. 80 g Normalized Microstresses, 0— B Direction, 

33 ’ 

9.72% Debonding, g e Loading 266 

yy 

A. 81 g Normalized Microstresses, 0— B Direction, 

11.11% Debonding, cr e Loading 267 

yy 

A.82 g Normalized Microstresses, 0— A Direction, 

12 ’ 

0.0% Debonding, <j e Loading 268 

yx 

A. 83 g Normalized Microstresses, O— A Direction, 

1.39% Debonding, cr e Loading 269 

yx 

A. 84 g Normalized Microstresses, 0— A Direction, 

12 ’ 

2.78% Debonding, cr e Loading 270 

yx 
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A.85 a Normalized Microstresses, 0— A Direction, 
12 

4.17% Debonding, a e Loading 
yx 

A. 86 a Normalized Microstresses, 0— A Direction, 

12 

5.56% Debonding, a e Loading 
yx 

A. 87 a Normalized Microstresses, 0— A Direction, 

12 

6.94% Debonding, a e Loading 
yx 

A. 88 a Normalized Microstresses, 0— A Direction, 

12 

8.33% Debonding, a e Loading 
yx 

A.89 <7 Normalized Microstresses, 0— A Direction, 

12 

9.72% Debonding, or® Loading 
yx 

A. 90 a Normalized Microstresses, 0— A Direction, 
12 

11.11% Debonding, <r e Loading 

yx 

A.91 a Normalized Microstresses, 0— B Direction, 
12 

0.0% Debonding, a e Loading 
yx 

A. 92 a Normalized Microstresses, 0— B Direction, 
12 

1.39% Debonding, <7 e Loading 
yx 

A.93 a Normalized Microstresses, 0— B Direction, 
12 

2.78% Debonding, a e Loading 
yx 

A.94 a Normalized Microstresses, 0— B Direction, 
12 

4.17% Debonding, a e Loading 
yx 
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279 
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XXVI 



A. 95 a Normalized Microstresses. 0— B Direction, 

5.56% Debonding, cr e Loading 

yx 

A.96 a Normalized Microstresses. 0-B Direction, 
12 5 

6.94% Debonding, <r e Loading 
y* 

A. 97 a Normalized Microstresses, O— B Direction, 

12 1 ’ 

8.33% Debonding, <j e Loading 

yx 

A. 9 8 a Normalized Microstresses, 0— B Direction, 

12 1 ’ 

9.72% Debonding, <r e Loading 

yx 

A.99 a Normalized Microstresses, 0— B Direction, 
12 

11.11% Debonding, cr e Loading 

yx 

A. 100 a Normalized Microstresses, 0— C Direction, 

12 

0.0% Debonding, cr e Loading 

yx 

A.101 a Normalized Microstresses, O-C Direction, 

IdS 

1.39% Debonding, <j e Loading 

yx 

A.102 a Normalized Microstresses, 0— C Direction, 

12 } 

2.78% Debonding, Loading 

yx 

A. 103 a Normalized Microstresses, 0— C Direction, 

12 1 5 

4.17% Debonding, a e Loading 

A.104 <j Normalized Microstresses, 0— C Direction, 

5.56% Debonding, a* Loading 

yx 
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A.105 

<j Normalized Microstresses, O-C Direction, 



6.94% Debonding, <r e Loading 

yx 
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A. 106 

a Normalized Microstresses, O-C Direction, 
12 



8.33% Debonding, cr e Loading 

yx 

292 

A. 107 

a Normalized Microstresses, O-C Direction, 
12 7 



9.72% Debonding, a e Loading 
y x 

293 

A.108 

a Normalized Microstresses, O-C Direction, 
12 7 



11.11% Debonding, a e Loading 

yx 

294 

A.109 

a Normalized Microstresses, 0— A Direction, 

23 7 



0.0% Debonding, <r e Loading 

yz 

295 

A. 110 

a Normalized Microstresses, O— A Direction, 

23 



1.39% Debonding, <r e Loading 

yz 

296 

A. Ill 

a Normalized Microstresses, O— A Direction, 

23 



2.78% Debonding, a e Loading 

yz 

297 

A. 112 

a Normalized Microstresses, O— A Direction, 

23 



4.17% Debonding, a e Loading 

yz 

298 

A.113 

a Normalized Microstresses, 0— A Direction, 

23 



5.56% Debonding, tr e Loading 
yz 

299 

A.114 

a Normalized Microstresses, 0— A Direction, 

23 



6.94% Debonding, cr e Loading 
yz 

300 
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A. 115 a Normalized Microstresses, 0— A Direction. 

23 

8.33% Debonding, <r e Loading 
yz 

A.116 <j Normalized Microstresses, 0— A Direction, 

23 

9.72% Debonding, <r e Loading 

yz 

A.117 a Normalized Microstresses, 0— A Direction, 

23 

11.11% Debonding, <j e Loading 

yz 

A. 118 <7 Normalized Microstresses, 0— B Direction, 

23 7 

0.0% Debonding, a e Loading 

yz 

A. 119 a Normalized Microstresses, 0— B Direction, 

23 

1.39% Debonding, <r e Loading 
yz 

A. 120 a Normalized Microstresses, O—B Direction, 

23 7 

2.78% Debonding, a e Loading 

yz 

A.121 a Normalized Microstresses, 0— B Direction, 

23 7 

4.17% Debonding, <7 e Loading 

yz 

A. 122 <j Normalized Microstresses, 0— B Direction, 

23 

5.56% Debonding, cr e Loading 

yz 

A. 123 a Normalized Microstresses, 0— B Direction, 

23 

6.94% Debonding, cr e Loading 

yz 

A. 124 a Normalized Microstresses, 0— B Direction, 

23 

8.33% Debonding, j e Loading 

yz 
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A. 125 a Normalized Microstresses. O—B Direction, 

23 

9.72% Debonding, a e Loading 

yz 

A. 126 c Normalized Microstresses, O—B Direction, 

23 

11.11% Debonding, <r e Loading 


A. 127 a Normalized Microstresses, 0— C Direction, 

23 

0.0% Debonding, <r e Loading 

yz 

A. 128 a Normalized Microstresses, O-C Direction, 

23 

1.39% Debonding, a e Loading 
yz 

A.129 a Normalized Microstresses, 0— C Direction, 

23 

2.78% Debonding, <7 e Loading 
yz 

A. 130 a Normalized Microstresses, 0— C Direction, 

23 

4.17% Debonding, cr e Loading 

yz 

A.131 a Normalized Microstresses, 0— C Direction, 

23 

5.56% Debonding, <r e Loading 

yz 

A.132 a Normalized Microstresses, 0— C Direction, 

23 

6.94% Debonding, <r e Loading 
yz 

A.133 a Normalized Microstresses, 0— C Direction, 

23 

8.33% Debonding, a* Loading 

yz 

A. 134 a Normalized Microstresses, O-C Direction, 

23 

9.72% Debonding, <r e Loading 

yz 
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312 


313 


314 


315 


316 


317 


318 


319 


320 


XXX 



A. 135 or Normalized Microstresses. O—C Direction, 

23 ’ 


11.11% Debonding, <r e Loading 321 

yz 

A. 136 a Normalized Microstresses, O—C Direction, 

0.0% Debonding, Thermal Loading 322 

A. 137 Normalized Microstresses, O-C Direction, 

I. 39% Debonding, Thermal Loading 323 

A. 138 a Normalized Microstresses, O—C Direction, 

2.78% Debonding, Thermal Loading 324 

A.139 a Normalized Microstresses, O—C Direction, 

4.17% Debonding, Thermal Loading 325 

A. 140 Normalized Microstresses, O-C Direction, 

5.56% Debonding, Thermal Loading 326 

A. 141 a ■ Normalized Microstresses, O—C Direction, 

6.94% Debonding, Thermal Loading 327 

A.142 a Normalized Microstresses, O—C Direction, 

8.33% Debonding, Thermal Loading 328 

A. 143 a Normalized Microstresses, O—C Direction, 

9.72% Debonding, Thermal Loading 329 

A.144 a Normalized Microstresses, O—C Direction, 

II. 11% Debonding, Thermal Loading 330 

A. 145 a Normalized Microstresses, 0— A Direction, 

0.0% Debonding, Thermal Loading 331 

A. 146 cr Normalized Microstresses, O— A Direction, 

1.39% Debonding, Thermal Loading 332 


xxxi 



A.147 

a Normalized Microstresses. 0— A Direction, 

31 

2.78% Debonding, Thermal Loading 

333 

A. 148 

a Normalized Microstresses, 0— A Direction, 

31 

4.17% Debonding, Thermal Loading 

334 

A.149 

a Normalized Microstresses, 0— A Direction, 

31 

5.56% Debonding, Thermal Loading 

335 

A. 150 

a Normalized Microstresses, 0— A Direction, 

31 

6.94% Debonding, Thermal Loading 

336 

A.151 

a Normalized Microstresses, 0— A Direction, 

31 

8.33% Debonding, Thermal Loading 

337 

A. 152 

cr Normalized Microstresses, 0— A Direction, 

31 

9.72% Debonding, Thermal Loading 

338 

A.153 

a Normalized Microstresses, 0— A Direction, 

31 

11.11% Debonding, Thermal Loading 

339 

A. 154 

<7 Normalized Microstresses, 0— C Direction, 

31 t 

0.0% Debonding, Thermal Loading 

340 

A.155 

a Normalized Microstresses, 0-C Direction, 

31 

1.39% Debonding, Thermal Loading 

341 

A. 156 

a Normalized Microstresses, 0— C Direction, 

31 

2.78% Debonding, Thermal Loading 

342 

A.157 

a Normalized Microstresses, 0— C Direction, 
31 

4.17% Debonding, Thermal Loading 

343 

A.158 

<7 Normalized Microstresses, 0-C Direction, 

31 

5.56% Debonding, Thermal Loading 

344 

A. 159 

a Normalized Microstresses, 0— C Direction, 

31 

6,94% Debonding, Thermal Loading 

345 
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A.160 
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A.163 
A. 164 
A. 165 
A. 166 
A.167 
A. 168 
A. 169 
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A. 171 


01 Normalized Microstresses, 0— C Direction, 
8.33% Debonding, Thermal Loading 

^ Normalized Microstresses, O-C Direction, 
9.72% Debonding, Thermal Loading 

<r Normalized Microstresses, O-C Direction, 
11.11% Debonding, Thermal Loading 

* Normalized Microstresses, O-C Direction, 
0.0% Debonding, Thermal Loading 

Normalized Microstresses, O-C Direction, 

I. 39% Debonding, Thermal Loading 

Normalized Microstresses, O-C Direction, 
2.78% Debonding, Thermal Loading 

a 12 Normalized Microstresses, O-C Direction, 
4.17% Debonding, Thermal Loading 

* Normalized Microstresses, O-C Direction, 
5.56% Debonding, Thermal Loading 

Normalized Microstresses, O-C Direction, 
6.94% Debonding, Thermal Loading 

Normalized Microstresses, 0— C Direction, 
8.33% Debonding, Thermal Loading 

<t Normalized Microstresses, O-C Direction, 
9.72% Debonding, Thermal Loading 

* Normalized Microstresses, O— C Direction, 

II. 11% Debonding, Thermal Loading 
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CHAPTER 1 
INTRODUCTION 

Composite materials have increased in importance as an engineering 
material for aerospace applications due to their high strength-to-weight ratio. 
Composite materials are increasingly being used in high temperature 
applications. To meet the need for high temperature composite materials, metal 
matrix and ceramic matrix materials have been utilized. Analysis and evaluation 
of these materials and their properties are very important if reliable designs are 
to be produced. 

In recent years the finite element procedure has been shown to be a 
viable method of structural analysis. A useful alternative in the evaluation of 
structures with repeated geometry is substructuring using the superelement 
method. Substructuring involves the partitioning of a finite element mesh into 
separate collections of elements called substructures. Each substructure is solved 
separately and then combined. Using superelements is virtually the same as 
substructuring, both in its mathematical implementation and in concept, with 
the primary exception being the user interface [1]. The superelement technique 
involves defining an image by copying, rotating, or mirroring a portion of a 
conventional finite element mesh ("conventional" shall imply a 
non-superelement approach). This image can be described as a structural 
building block and may be repeatedly included as part of the overall structure. 
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No modeling data is required for the description of the image other than a list of 
its boundary points and information on its orientation with respect to the initial 
representation [2]. 

Traditionally, fibrous composite materials have been modeled as 
structures with evenly spaced fibers surrounded by a matrix material, the fibers 
having consistent diameters and running parallel throughout the composite. 
Though, in reality, some deviations from this geometry may be present, past 
work has shown that the assumption of consistent spacing is a good one [3] [4]. 

The repeated geometry of the idealized composite readily lends itself to 
representation by the superelement technique. By conventionally modeling a 
unit cell of composite material and then imaging that unit cell, a composite 
structure of desired dimensions may be described. Analysis by the superelement 
technique is quite efficient and more easily implemented than conventional finite 
element methods. 

The composite superelement mesh is of modular nature, consisting of a 
unit cell (figure 1.1) and images of the unit cell; therefore, discontinuities can 
easily be added by substituting a conventional mesh, which models the 
discontinuity, in place of one of the images. This conventional mesh can also be 
imaged to simulate multiple discontinuities within the composite material. Some 
discontinuities of interest that can be modeled in this manner are partial bonding 
and fiber fracture, both common in metal matrix composite materials. The 
modular nature of the superelement model also allows for efficient representation 
of hybrid and aligned short fiber composite materials. The hybrid composite can 
be modeled by changing the fiber material properties of two separate primary 
superelements and using these primaries and their images to build the desired 
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structure. An aligned short fiber composite may be formed by including sections 
of pure matrix between lengths of fibers. The generation of a bead— filled 
composite may be accomplished using similar methods. 

Once a superelement mesh is determined to be void of modeling 
inaccuracies, it can then be altered to investigate new areas of composite 
technology without the high cost of experimental analysis. Some of these areas 
may be the use of hollow reinforcing fibers or fibers with non— circular cross 
section. The effect of variations in fiber diameter along the fiber’s length could 
also be investigated. 
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CHAPTER 2 
OBJECTIVE 


The objective of this paper is to show the degradation of various 
properties of metal matrix composite materials in relation to the amount of fiber 
disbond. The work also shows the mechanics of load transfer from the 
surrounding material to the debonded fiber. Since the debonding is only for a 
fraction of the length of the fiber, the fiber becomes structurally active in the 
bonded region and a gradient of load transfer from the matrix to the fiber can be 
witnessed. In the region where there is no connectivity between the fiber and the 
surrounding matrix, it is expected that the fiber will deform rigidly relative to 
the rest of the structure; therefore, theoretically, there will be no stresses in this 
region of the fiber. 
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CHAPTER 3 


BACKGROUND 

The composite system considered for this work is 
P100— Graphite/ Copper. Analyses for high temperature and room temperature 
constituent material properties were considered. The utilization of copper for the 
metal matrix material is significant because of its relatively low thermal 
resistivity coefficient, which allows for a high rate of heat transfer. Thermal 
stresses in metal matrix composites can be very high due to the difference in the 
thermal expansion coefficients of the constituents. The high thermal 
conductivity of the copper matrix allows for the dissipation of heat energy, 
causing a decrease in thermal straining and a subsequent lowering of the stress 
state in the composite. 

A problem arises when choosing properties to describe each of the 
constituent materials because the state of the material, and thus its properties, 
depends on the history of the material. An important factor in the material’s 
behavior is its rate of cooling during the fabrication process. The rate of cooling 
determines the degree of crystallization. If copper is quickly cooled it will tend to 
be more amorphic, with many small crystals, and therefore more ductile. If it is 
cooled slowly it will have time for larger crystals to form and will be more brittle 
in nature [5]. It is difficult to test the in situ state of the copper after fabrication 
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of the composite; therefore, average material properties for the copper matrix 
must be assumed. This assumption may be a source of error when comparing 
finite element results to actual test data. 

The superelement method is a form of substructuring used in the 
MacNeal Schwendler Corporation (MSC) version of NASTRAN where the 
computer takes on most of the bookkeeping burden as well as doing the entire 
analysis in one computer run. For this study Version 61 of MSC NASTRAN was 
used in the superelement analysis. All NASTRAN finite element work was done 
on a Cray-XMP supercomputer utilizing a solid state storage device. The 
database for each run was saved on the Cray and restarted using the next set of 
boundary conditions. This method allowed for faster computer turnaround 
during analysis. 

The cases chosen are but a few of a multitude of possible partial 
bonding configurations. However, these cases are sufficient to demonstrate the 
load transfer and stress concentration effects of this phenomena and help to 
understand the effects of partial bonding on the composite’s material properties. 
In the cases examined here, all fibers are running parallel to each other with their 
longitudinal axes in the X direction. Methods exist for transforming the material 
properties calculated from this fiber orientation to other material axis systems 
[«]• 

The mesh that was employed in this work consisted of a cluster of nine 
fibers in a three by three matrix (figure 3.1) and represents a fiber volume ratio 
of 0.466. The initial version of the mesh was fabricated with the center cell 
modeled as a primary superelement and all surrounding cells were images of that 
primary. This mesh was utilized when no fibers were debonded and when all nine 
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fibers were debonded equally. Another mesh was formulated for when only the 
center fiber of the nine cell model was debonded. For this case the center cell 
was modeled using conventional finite element methods. The primary 
superelement was defined to be to one side of the center fiber and the remaining 
perimeter superelements were modeled as images of this one. Minor 
modifications were made to each of the meshes to allow for the modeling of 
varied amounts of disbond between fiber and matrix. 

The methodology utilized to include the disbond is described briefly 
here. Duplicate grid points were defined for all grids around the circumference of 
the fiber, at the interface between fiber and matrix. This duplicate grid point 
was included in the analysis only when the portion of fiber where it resided was 
to be debonded from the matrix. The debonding took place when the 
connectivity cards that define the perimeter of the fiber were altered so that one 
of the duplicate grids was associated with the fiber, and one with the matrix. No 
connectivity then exists across the interface, producing a crack of zero width, 
with fiber on one side and matrix on the other (figure 3.2). In the work described 
here, the total circumference of a fiber was released together. Future work may 
include analyses where the fiber is not totally circumferentially disbonded, but 
this was not approached here. The amount of disbond was varied by debonding 
different lengths of fiber as described above. No attempt was made to find the 
loading at which the debonding will take place. Instead, interest was placed on 
the effect of debonding on the composite’s material properties. 
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Figure 3.2 — Modeling of Fiber Disbond 
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In the model adopted for this work, thin elements were placed around 
the circumference of the fibers (figure 1.1). For the work discussed here, this thin 
layer was defined to be matrix material but was included for future work 
intended for this mesh. These elements have large aspect ratios but have been 
shown not to significantly effect the material property calculations under the 
loading conditions used. This layer allows for the ability to include an interphase 
region which is common in metal matrix composites. This interphase region is 
formed in situ from a chemical reaction between the fiber material and matrix 
material during the fabrication process [7]. It can also be used to model a coating 
on the fiber which can decrease areas of high stress. Reference [8] provides a 
discussion on this effect. This interface is commonly the most highly stressed 
area of a composite material. It is advantageous to lower the stress 
concentrations in these areas by placing a material of intermediate modulus or a 
ductile one between the fiber and the matrix. When a graded modulus interphase 
is employed, the interphase layer is defined to be one of intermediate modulus 
and the modulus ratio of any two neighboring components is lowered thus 
lowering stress concentrations due to modulus mismatch. If the interface region 
is defined to be a ductile material then local deformation capability is built into 
this area of stress risers resulting in at least partial dampening of the stress 
concentration effects due to modulus mismatch. 

The following are the procedures used to calculate ply material 
properties from the finite element output. These methods are consistent with the 
mechanics of materials approach to solid structures [9]. 
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u , v , and E are obtained from the same loading conditions. In 
this case the front face is fixed in the X (u^ = 0.0) and the back face is displaced 
in the X (u g = u) (figure 3.3). Lines of geometric symmetry on the top and side 
faces are restricted from motion in such a way to insure symmetry once the 
model is loaded. 

E is calculated by finding the total force over the displaced face and 
dividing this force by the area of this face yielding an equivalent applied stress 
(figure 3.4). 


a 6 

in 


E F 



3.1 


The strain is then calculated by dividing the applied displacement by the length 
of the specimen. 


a & is divided by e to equal Young’s modulus in the 11 direction, 
ill ill ° 


E 

in 



3.3 


v is found by finding the average deflection in the Y direction as a 
result of the enforced displacement in the X. Dividing this deflection by the 
width of the model yields a strain in the 22. 
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Figure 3.4 — Model Dimensions 
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The Poisson’s ratio, v , 

112 


is calculated from equation 3.5. 
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v is found by similar methods. 
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v and E are both calculated from the same loading conditions. 
The face with the negative Y direction as its normal is fixed in the Y direction 
(v = 0.0) (figure 3.3) and the face with the positive Y direction as its normal 

3 * 

has an enforced displacement in the positive Y direction (v = v) resulting in 

g 

tension transverse to the fiber direction. The equivalent applied stress is then 
calculated by averaging the resultant forces over the face and dividing by the 
face’s area. 
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In calculating E the strain in the 22 is required. The strain is 
calculated by dividing the applied displacement by the width of the specimen. 


The equivalent stress divided by this strain is equal to Young’s modulus in the 22 
direction. 



122 
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v is found by finding the average deflection in the Z direction as a 

1 2 3 

result of the enforced displacement in the Y. Dividing this deflection by the 
height of the model yields the strain in the transverse direction (33). 
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The ratio of the strain in the 33 over the strain in the 22 yields the Poisson’s 
ratio in the 23. 
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E 33 and v are calculated from another unique set of boundary . 
conditions. In this case, the face whose normal is in the negative Z direction is 
fixed in the Z (w^ — 0.0) (figure 3.3), and the face whose normal is in the positive 
Z direction is displaced in the Z (w = w). The equivalent applied stress is 
computed by averaging the forces on the face where the deflection was applied 
and dividing by the area of that side. 


E F 

<7 e = T — 

133 A 

3 


3.13 


E l 33 * S caJcu * atec * by dividing this equivalent stress by the strain in the 
33 direction. The strain in the 33 is ascertained by dividing the applied 
displacement in the Z by the height of the specimen. 


The equivalent stress is divided by this strain and is equal to Young’s modulus in 
the 33 direction. 


E = 

133 


& c 
133 

£ 

133 


3.15 


The first step in finding u is obtaining the average deflection in the Y 
direction as a result of the enforced displacement in the Z. Dividing this 
deflection by the width of the model yields a strain in the 22. 
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3.16 


6 \22 


The ratio of this strain over the strain in the 33 yields the Poisson’s ratio, v 

132 


—e 

v = 111 . 

133 


123 


3.17 


For the determination of G , an enforced displacement in the X 
direction was placed on the face with the positive Y axis as its normal. The 
shear strain was then calculated by dividing this displacement by the width of 
the model. 


6 = 
121 


U 


3.18 


a 


A FORTRAN program was used to glean and sum the resulting grid forces on 
this face from the NASTRAN output deck to achieve an equivalent applied force. 
The effective shear stress in the 21 direction is then calculated by dividing this 
force by the area over which it is applied. 


E F 


<7 e = 

12.1 


121 


3.19 


The shear strain for small displacements is given by the enforced deflection 
divided by the width of the specimen. G is obtained by dividing the shear 
stress in the 21 direction by the shear strain. 
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3.20 


G 

121 


a e 
121 
e 

121 


G is found by similar methods. For this calculation the face with the 
positive Z axis for its normal is used for the location of the enforced 
displacement. 


c i3i s 

3 


E F 

= r 131 

131 A 
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3.21 


3.22 


G 

131 


0 * 

131 

€ 

131 


3.23 


G i 23 can be obtained b y applying an enforced displacement in the 
positive Z direction to the side whose normal is the positive Y axis. The total 
resulting force on this face is calculated from the finite element output. This 
total shear force in the 23 direction is divided by the area of the side over, which 
it was applied, resulting in the average shear stress in the 23 direction. 

E F 

«I,, = TT* 3.24 

2 
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The small displacement shear strain in the 23 is then calculated by dividing the 
applied displacement in the Z direction by the width of the model. 



3.25 


G is calculated by dividing the effective shear stress, o* , by e . 
123 123 123 



3.26 


ol , a , and a are all found by the same loading conditions. In 

111 ’ 122 ’ 133 

this case a plane in the center of the model, with the positive X direction acting 

as its normal, is restricted from movement in the X direction (u = 0.0). This 

restriction assures symmetry about the center of the model. Again, lines of 

symmetry are enforced on the sides of the model. A thermal load is then applied 

to the model, T = T . 

0 

The first step in finding o is determining the average of the 
displacements on the end with the positive X axis as its normal and dividing by 
the length of the model, s , resulting in the strain in the 11 due to the change in 
temperature. 



3.27 
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Now, by dividing the strain by the change in temperature, one can obtain the 
thermal expansion coefficient in the 11 direction. 


a 

in 



3.28 


& l22 is found by taking the average displacement in the Y direction and 
dividing by the change in temperature and the width. 

V = "(VT)(s 2 ) 329 

The calculation of a is similar to that of a , with the exception 

133 122 

being that the average displacement in the Z is divided by the change in 
temperature and by the height of the model [10]. 

“i2 2 = iVTflsj 330 

3 

In all of the previous loading conditions boundary conditions were 
applied to enforce displacement symmetry. On each face of the above loading 
conditions and for each degree of debonding investigated, normalized constituent 
microstresses in three directions, each starting in the center of the model and 
moving toward the perimeter, were determined. These lines can be referenced 
easily by referring to the line segments defined in each of the three directions. 
The first line is from the middle of the center fiber to the edge of the mesh in the 
positive Z direction (line segment O-A). The second starts in the center of the 
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center fiber and diagonally bisects the angle between the positive Y and positive 
Z directions, forming a 45 degree angle with each and lying in the same plane 
(line segment O-C). The last line of stresses, again, starts at the middle point of 
the centermost fiber and ends at the edge of the mesh running in the positive Y 
direction (line segment 0— B) (figure 3.5). A FORTRAN code was written to 
acquire constituent microstresses in the 11, 22, 33, 12, 23, and 31 directions along 
these lines and normalize them with respect to the equivalent applied stress. The 
normalized microstresses were gathered into tables and, if they were found to be 
significant relative to the ultimate stress of the constituent material in which 
they appeared, they were plotted. These plots allow for the investigation of areas 
of high stress due to debonding. These microstresses are not obtainable 
experimentally because they are internal to the structure. 
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Figure 3.5 — Locations of Stress Plots 



CHAPTER 4 
THEORETICAL 


The superelement method is a form of substructuring used in MSC 
NASTRAN. The important difference between substructuring and the use of 
superelements is that in the use of superelements most of the logistics of 
maintaining the substructure is done by the program and not the user. 
Mathematically, the use of superelements is equivalent to substructuring [11, 12, 
& 13]. A brief description of substructuring theory follows. 

Consider the static analysis equation 

w (",)-{ p ,i u 

where £ j is the stiffness matrix, | | is displacement vector and the vector 

| P f j is the load vector. The "f" subscript indicates statically independent or 
"free" variables and reflects the elimination of degrees of freedom due to 
multipoint and single— point constraints. These matrices can be partitioned into 
the a-set and the o-set where the o-set includes all interior degrees of freedom of 
the superelement and the a-set includes all exterior degrees of freedom. This is 
given by the equation 
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i 

i 

I 


kiki (".} i p .i 

4.2 

[ K ..][ E «]j [(".)) |f p .t 

Multiplying yields 

[ E Jkl + [ K JI“.)-l p .) « 


K]l".) + k.]KM p .) 

The bar over a variable indicates the quantity is a partition of its parent matrix. 
Solving the second equation for the o-set variables yields 

{■.}- [ K ~]'{ p .}-[ K »] , [ K „]{".} « 

Defining new terms we have 

where 



25 



kl'Kl'k) 


[ G J-K] [ K J 

The quantity, | u° J, is the partial solution obtained when loads are placed on 
interior points and external points are constrained. If no loads exist on points 
interior to the superelement then j u° J is null. Every superelement in the 
structure has a j u° j vector. The other part of the interior displacement 
solution is found by combining equations 4.8 and 4.4 yielding 

kkKlk) 49 

which is the displacement due to boundary point motion. Again, a j u* J vector 
is found for each superelement in the structure. The total solution is found by 
superposition of these two partial solutions. Combining the above equations 
yields equation 4.10. 

[ K Jk)-{ p .) 


where 
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The J J matrix is the stiffness of the superelement with respect to its 
boundary points and exists for each superelement. The j P j vector is 
representative of the boundary loads that are transferred to the rest of the 
structure. 

These partial solutions are stored and later added to the remainder of 
the structure. It is these matrices that, under conventional substructuring, 
would have to be manipulated by the user, whereas in the superelement 
technique the computer maintains this database. 

The mechanics of material approach to calculating material properties 
from finite element output is based on the definitions of the properties 
themselves. Boundary conditions are applied so as to do away with any ill effects 
of the loading that might change the calculated property. Bending is one such ill 
effect. A description of the utilization of the mechanics of material approach for 
the calculation of composite material properties was addressed earlier. 


CHAPTER 5 


PROCEDURE 

The initial investigation into partial bonding was preceded by some test 
cases to confirm the accuracy of the methods used to obtain the material 
properties from the finite element output. To do this, the mesh was executed 
using a monolithic material, P100— graphite (the fiber material). The loading 
cases discussed earlier were applied to the model and the monolithic material 
properties were back-calculated. These properties were compared to those of 
P100— graphite used as input to the finite element code. The calculated material 
properties obtained from the finite element output were very close to those used 
as input (table 5.1). This confirms the viability of the this method of calculating 
composite material properties. 

A set of reference finite element runs were then executed. These runs 
contained no debonding. Material properties were calculated from the output of 
these finite element runs (table 5.2). Caruso and Chamis [14] calculated 
composite material properties using similar methods and compared them to those 
predicted by Hopkins and Chamis [15] with very good correlation. 

Debonding of a single fiber in the center of the nine cell model was now 
considered. Room temperature (70° F) constitutive material properties were 
used (table 5.3). In all loading conditions described here, fibers that were to be 
debonded did not have an applied displacement assigned to them. This allowed 
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CALCULATED MATERIAL PROPERTIES 
FOR ALL PlOO- GRAPHITE MODEL 


CONSTITUTIVE CALCULATED 

PlOO- GRAPHITE PlOO- GRAPHITE 

UNITS PROPERTIES PROPERTIES 


E 

111 

psi 

105.0 x 10 6 

105.0039201 x 10 8 

E 

122 

psi 

0.90 x 10 6 

0.900165249 x 10 8 

E 

133 

psi 

0.90 x 10 6 

0.900165249 x 10 6 

G 

1 12 

psi 

1.10 x 10 8 

1.100001462 x 10 6 

G 

123 

psi 

0.70 x 10 6 

0.6999990253 x 10 8 

G 

113 

psi 

1.10 x 10 8 

1.100001462 x 10 6 

V 

1 12 

in/ in 

0.200 

0.2000547581 

V 

123 

. 

in /in 

0.250 

0.2496567996 

V 

1 1 3 

in/ in 

0.200 

0.2000547581 

a 

1 l 1 

in/in/°F 

-0.90 x 10~ 8 

-0.90 x 10* 8 

a 

122 

in/in/°F 

5.60 x 10* 6 

5.60 x 10- 8 

a 

133 

in/in/°F 

5.60 x 10* 8 

5.60 x 10* 8 


Table 5.1 - Calculated Material Properties for all PlOO-Graphite Model 
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COMPOSITE MATERIAL PROPERTIES 
NO DEBONDING - ALL FIBERS LOADED 
PlOO- GRAPHITE / COPPER 


COMPOSITE 
UNITS PROPERTIES 


E 

l i i 

psi 

58.03x10* 

E 

122 

psi 

7.19x10* 

E 

133 

psi 

7.19x10* 

G 

112 

psi 

3.36x10* 

G 

123 

psi 

2.66x10* 

G 

1 1 3 

psi 

■ 

3.36x10® 

V 

1 12 

in/ in 

0.291 

V 

123 

in/ in 

0.246 

V 

1 13 

in/ in 

0.291 

a 

1 1 1 

in/in/°F 

j 1 . llxlO - * 

a 

122 

in/in/°F 

11.41x10-* 

Q 

133 

in/in/°F 

mm 


Table 5.2 - Composite Material Properties - No Debonding 
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ROOM TEMPERATURE CQNSTITUITIYE MATERIAL PROPERTIES 
PlOO- GRAPHITE / COPPER 

OBITS FIBER MATRIX 


E 

1 i 

psi 

105.0 x 10 6 

17.0 x 10 8 

E 

22 

psi 

0.90 x 10* 

17.0 x 10 8 

B 

33 

psi 

0.90 x 10 s 

17.0 x 10 8 

G 

12 

psi 

1.10 X 10‘ 

6.54 x 10‘ 

G 

23 

psi 

0.70 x 10 8 

6.54 x 10 8 

G 

1 3 

psi 

1.10 x 10 8 

6.54 x 10 6 

V 

12 

in/ in 

0.200 

0.300 

V 

23 

in/ in 

0.250 

0.300 

V 

1 3 

in/ in 

. 

0.200 

0.300 

a 

1 1 

in/in/°F 

-0.90 x 10‘ 8 

9.80 x 10' 6 

a 

22 

in/in/°F 

5.60 x 10‘ 8 

. 9.80 x 10’ 8 

a 

33 

in/in/°F 

5.60 x lO’ 6 

9.80 x 10" 6 


Table 5.3 - Room Temperature Constitutive Material Properties 
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these fibers to behave as if they had been fractured. Debonding dije to stress 
increases resulting from a broken fiber is a common occurrence in metal matrix 
composites. In this first series of runs, debonding was propagated from only one 
end of the model, producing an asymmetric disbond. Each element around the 
circumference of the fiber was allowed to debond using the methods described 
earlier. The debonding was done for each successive layer of elements along the 
length of the fiber. After the mesh was changed to reflect the next layer of 
debonding, the loading conditions were applied and material properties 
calculated. Results were plotted for each of the material properties versus 
percent disbond of the total mesh (figures 5.1 — 5.4). Least— squares linear 
regression was used to obtain equations for these lines. This allowed the results 
to be used as aid in the prediction of composite property degradation due to fiber 
debonding [16]. These plots show the properties affected most by debonding and 
to what degree they are affected. These results can also be used indicate partial 
bonding in cases where material properties deviated from predicted values. 
When suspecting the presence of partial bonding, closer attention should be paid 
to those properties found to be more sensitive to the presence of the fiber 
debonding. 

The next series of runs had the same geometry as those mentioned 
above except that high temperature constituent material properties were used 
along with the symmetric debonding (table 5.4). The material properties found 
in table 5.4 were obtained from the METCAN (Metal Matrix Composite 
Analyzer) computer code [17]. The temperature selected was 1500° F. Here, all 
debonding was done symmetrically with respect to the center of the model. The 
results of the material property calculations are shown in figures 5.5 through 5.8. 
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MODULUS - psi (xIO 6 ) 


EFFECT OF FIBER DEBONDING ON 
MODULUS - E1I1, E122, E133 



PERCENT DISBOND 


Figure 5.1 — Effect of Center Fiber Debonding on Modulus at Room Temperature 




POISSONS RATIO 


EFFECT OF FIBER DEBONDING ON 
POISSONS RATIO ~ i412, yll3, vffl, A32 



PERCENT DISBOND 

Figure 5 2 - Effect of Center Fiber Debonding on Poisson’s Ratio at Room Temperature 





SHEAR MODULUS - psi (x10 e ) 


EFFECT OF FIBER DEBONDING ON 
SHEAR MODULUS - GI21, G131, G123 



PERCENT DISBOND 

Figure 5.3 - Effect of Center Fiber Debonding on Shear Modulus at Room Temperature 




TEC (xicr 6 in/in/°F) 


EFFECT OF FIBER DEBONDING ON 
THERMAL EXPANSION COEFFICIENT (TEC) - ortll, a!22, a!33 



PERCENT DISBOND 

Figure 5.4 - Effect of Center Fiber Debonding on Thermal Expansion Coefficient at Room Temperature 




HIGH TEMPERATURE (1500° P) CONSTITUITIVE MATERIAL PROPERTIES 

P100- GRAPHITE / COPPER 

' 4 


UNITS TIBER MATRIX 


E 

1 1 

psi 

93.8 x 10“ 

8.87 x 10® 

E 

22 

psi 

0.904 x 10 s 

mm 

E 

33 

psi 

0.904 x 10 6 

8.87 x 10‘ 

6 

12 

psi 

0.982 x 10 6 

3.41 x 10 6 

G 

23 

psi 

0.625 x 10* 

3.41 x 10 6 

G 

13 

psi 

0.982 x 10® 

3.41 x 10 6 

V 

12 

in/ in 

0.179 

0.300 

V 

23 

in/in 

0.223 

0.300 

V 

1 3 

in /in 

0.179 

0.300 

a 

1 1 

in/in/°E 

mm 

19.55 x 10* 6 

a 

22 

in/in/°F 

6.27 x 10* 6 

19.55 x 10- 6 

a 

33 

in/in/°F 

6.27 x lO -6 

19.55 x 10’ 6 


Table 5.4 — High Temperature Constitutive Material Properties 
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MODULUS - psi (xIO 6 ) 


EFFECT OF FIBER DEBONDING ON 
MODULUS - Bill, E122, E133 
HIGH TEMPERATURE (1500 °F> 



PERCENT DISBOND 

Figure 5.5 - Effect of Center Fiber Debonding on Modulus at High Temperature 




POISSONS RATIO 


EFFECT OF FIBER DEBONDING ON 
POISSONS RATIO - i412, t/113, 1/123, vY 32 
HIGH TEMPERATURE (1500 °F) 



PERCENT DISBOND 


Figure 5.6 — Effect of Center Fiber Debonding on Poisson’s Ratio at High Temperature 




EFFECT OF FIBER DEBONDING ON 
SHEAR MODULUS - G121, G131, G123 
HIGH TEMPERATURE (1500 °F) 
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G131 

G123 


Figure 5.7 - Effect of Center Fiber Debonding on Shear Modulus at High Temperature 




TEC (x10' e ln/In/°F) 


EFFECT OF FIBER DEBONDING ON 
THERMAL EXPANSION COEFFICIENT (TEC) - cclll, 0(122, al33 
HIGH TEMPERATURE (1500 °F) 



PERCENT DISBOND 


Figure 5.8 - Effect of Center Fiber Debonding on Thermal Expansion Coefficient at High Temperature 





Note that these plots reflect a little over eleven percent debonding of the entire 
composite. This is a result of releasing only the center fiber for these cases 
considered here. The 11% debonding value shown on these plots represents the 
point at which the center fiber is totally debonded and the surrounding eight 
fibers are fully bonded to the matrix. This yields just over eleven percent (or one 
ninth) of the total circumferential area of all fibers not bonded to the matrix. 
When the center fiber is considered completely debonded it is still connected to 
the matrix material by a ring of nodes in the plane of symmetry. This residual 
attachment of the fiber’s circumference to the matrix material can add some 
stiffness transversely and in shear. 

The next series of conditions considered were when all nine fibers in the 
model were debonded symmetrically and room temperature conditions were 
imposed. Debonding was obtained by simultaneously removing the connectivity 
between fiber and matrix for two individual finite element lengths, one from each 
end of the model, to assure symmetry. In this case the percent debonding can go 
up to 100% of the total fiber length. Once again, the connectivity of a ring of 
nodes in the plane of symmetry around each fiber was maintained at 100% 
debonding. Material properties were calculated by applying each of the loading 
conditions described earlier. These results were plotted to show the effect of fiber 
debonding on each property when all fibers were equally debonded (figures 5.9 — 
5.12). As before, the least— squares method was used to fit an equation to each of 
these lines. 

Fiber debonding was then evaluated for high temperature conditions. 
All nine fibers were released as described above using material properties that 
reflected a use temperature of 1500° F. Material properties were calculated and 
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MODULUS - psi (xIO 6 ) 


EFFECT OF FIBER DEBONDING ON 
MODULUS - Elll, EI22, E133 



PERCENT DISBOND 

Figure 5.9 — Effect of Total Fiber Debonding on Modulus at Room Temperature 




POISSONS RATIO 


EFFECT OF FIBER DEBONDING ON 
POISSONS RATIO - t412, i/113, i423, i432 



PERCENT DISBOND 

Figure 5.10 - Effect of Total Fiber Debonding on Poisson’s Ratio at Room Temperature 






EFFECT OF FIBER DEBONDING ON 
SHEAR MODULUS - G121, G131, G123 
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Figure 5.11 - Effect of Total Fiber Debonding on Shear Modulus at Room Temperature 




TEC (xIO " 6 in/in/°F 


EFFECT OF FIBER DEBONDING ON 
THERMAL EXPANSION COEFFICIENT (TEC) - alll. al22, a!33 



PERCENT DISBOND 


Figure 5.12 - Effect of Total Fiber Debonding on Thermal Expansion Coefficient at Room Temperature 





plotted for varying degrees of fiber debonding ranging from 0% to 100% 
debonding (figures 5.13 — 5.16). These results show the effect of high 
temperature on composite material properties with fibers debonding in groups. 
Least— square equations describing the material property changes due to this 
effect were calculated and are discussed later. 



MODULUS - psi (xIO 6 ) 


EFFECT OF FIBER DEBONDING ON 
MODULUS - Bill, E122, E133 
HIGH TEMPERATURE (1500 °F) 



PERCENT DISBOND 

Figure 5.13 — Effect of Total Fiber Debonding on Modulus at High Temperature 



POISSONS RATIO 


EFFECT OF FIBER DEBONDING ON 
POISSONS RATIO - i/U2, 1/113, 1/123, 1/132 



PERCENT DISBOND 

Figure 5.14 - Effect of Total Fiber Debonding on Poisson’s Ratio at High Temperature 





SHEAR MODULUS - psi (xIO 6 ) 


EFFECT OF FIBER DEBONDING ON 
SHEAR MODULUS - G121, G131, G123 
HIGH TEMPERATURE (1500 °F) 



PERCENT DISBOND 


Figure 5.15 - Effect of Total Fiber Debonding on Shear Modulus at High Temperature 
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EFFECT OF FIBER DEBONDING ON 
THERMAL EXPANSION COEFFICIENT (TEC) - aUl, al22, al33 
HIGH TEMPERATURE (1500 ®F) 



PERCENT DISBOND 

Figure 5.16 — Effect of Total Fiber Debonding on Thermal Expansion Coefficient at High Temperature 




CHAPTER 6 


RESULTS 

The degradation in the following material properties is influenced by the 
difference in the corresponding fiber and matrix material properties. As the 
fibers become structurally inactive then the matrix material properties will have 
a more significant role in determining the composite’s global properties. If the 
difference between fiber and matrix properties is not great then little change will 
be seen globally. Therefore there are two parts in assessing the sensitivity of a 
composite’s material properties to fiber debonding: 1) the physical interaction 
between fiber and matrix as stresses are transferred between the two and 2) the 
relative different between the constituent’s material properties. Therefore, while 
viewing the following results, consideration must be given to the fact that 
although a property is declared to be insensitive to debonding using these 
materials (P 100/Copper) the apparent insensitivity may be due to a small 
difference in constituent material properties. In some cases, it will be seen that a 
composite material property can be lower than that of the lowest constituent 
property. 
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6.1) Room Temperature Moduli with Center Fiber Debonding 

Consider the results portrayed in figure 5.1. This figure shows the 
degradation of the modulus in the 11, 22 and 33 directions due to debonding of 
the center fiber only. These results reveal that the longitudinal modulus has 
decreased by about 8% while debonding was at 11.1% of the total circumferential 
fiber area in the model. When the least squares function approximation is used 
to quantify this line the resulting equation is found to be 

= 5.706*10 7 - (4.368x10 5 )P 6.1 

where "P M is the percent of circumferential fiber area that is debonded. 

E i22 and E were solved for by utilizing the results from different loading 
conditions but were consistently calculated to be virtually equally to each other 
as shown by the coincident lines on figure 5.1. Because of the geometry of the 
idealized composite, properties in the 22 and 33 directions should be equivalent 
[17]. The coincidence of these lines suggest the validity of this approach. These 
lines can be approximated by the equation 

E = E = 7.191x10® - (9.499 x10 3 )P 6.2 

lwu 

The percent decrease in E^ and E for 11.1% total fiber area debonding is 

1.5%. From these results it can be concluded that E is more sensitive to fiber 

in 
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debonding than the transverse moduli. Even stilly the decrease ( in the 

longitudinal modulus is reasonably small considering E is 6 times greater than 

E . From further results it will become evident that E is the property that is 
m m 

the second most sensitive to fiber debonding. 

6.2) Room Temperature Poisson’s Ratio with Center Fiber Debonding 

The effect of fiber debonding on Poisson’s ratio for center fiber debonding 
is depicted in figure 5.2. Once again there is consistency with expected results 
since v coincides with v and v coincides with v and each were solved 

112 113 123 132 

using the output from different loading conditions. The least squares 
approximation for v and u is given by 

v = v = 0.293 - (3.590*10" 4 )P 6.3 

112 113 v 

When the center fiber is completely debonded (11.1% total fiber length 
debonding), v and v have degraded by 0.5%. This is a very small amount 

112 113 

and may be considered insignificant implying that v and v are quite 

112 113 

insensitive to fiber debonding. This insensitivity may be due to the small 

difference between v and v ( v is 1.5 times larger than v ). 

f 12 m m fl2 

v and v have a least squares approximation of 

123 132 


!/=!/= 0.246 - (3.327x10' 4 )P 6.4 

123 132 
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Since v is 1.2 times greater than v and a 1.6% decrease in v is witnessed 

m f23 123 

when total center fiber debonding is achieved; it can be concluded that geometric 
effects dominate over the effect of the relative increase in the role of the higher 
matrix Poisson’s ratio. From these results we can truly say that v is more 

123 

sensitive to fiber debonding than i/ . Given the same amount of debonding, 

there was relatively less degradation of v than v even though the difference 

112 123 

between the constituent’s respective Poisson’s ratios was greater. 

6.3) Room Temperature Shear Moduli with Center Fiber Debonding 

The decrease in shear moduli due to center fiber debonding is shown in 
figure 5.3. As can be expected, G and G are coincident, even though they 
are calculated from the results of different loading conditions. G decreased by 
3% as a result of complete debonding of the center fiber. Note that G is six 

m 

times larger than G^ and yet there is a decrease in G as increasing amounts of 
the center fiber become structurally inactive. The least squares linear 
approximation of the G^ and G is given by 

G = G = 3.357*10 6 -(9.495xl0 3 )P 6.5 

112 113 V ' 


! 

i 

j 


Figure 5.3 also shows the effect of center fiber debonding on G . The 

123 

value of G decreased by 2% when the center fiber was completely debonded, a 
decrease of structurally active fiber by 11.1%. Note again that though G is 

f 23 

nine times less than G , there is a decrease in G as more of the fiber is 

m 123 
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debonded. Matrix properties would be expected to play a more significant role in 
the determination of the composite’s global properties as increasing amounts of 
the fiber become structurally inactive. The approximate equation for the 
degradation of G is given by 

123 

G = 2.720*10 6 - (4.780*10 3 )P 6.6 

123 

G , G as well as G can be considered to be relatively insensitive to 

112 * 113 123 

the debonding of the center fiber due to the small changes in these material 
properties. 

6.4) Room Temperature Thermal Expansion Coefficients with Center Fiber 
Debonding 

The thermal expansion coefficients are plotted in figure 5.4. The lines 
representing the transverse thermal expansion coefficients are coincident, as to be 
expected. The debonding conditions considered here produced a decrease in the 
transverse thermal expansion coefficients by 0.3%. This decrease is minimal and 
may be considered to be insignificant. Note also that a is only 1.75 times 

a 

greater than a . The least squares approximation of and is given by 

a = a = 1.141*10- 5 +(3.003 x10* 9 )P 6.7 

122 133 
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The line describing the change in the longitudinal thermal expansion 
coefficient is also displayed in figure 5.4. increased by 17.1% when the center 
fiber was totally debonded, which is equivalent to 11.1% debonding of the total 
fiber area. The longitudinal thermal expansion coefficient for the fiber material 
is negative but a is positive. This difference may be the cause of such a 

m 

dramatic increase in a . Of all the material properties considered so far. a 

in in 

seems to be the most sensitive to center fiber debonding. This increase in 
may be considered significant. The least squares linear approximation of the line 
depicting the increase in on figure 5.4 is given by 

a = 1.226xl0” 6 + (1.362*10" 8 )P 6.8 

Analysis of the curve depicting the increase in indicates that most of the 
increase occurs initially and then the line tends to decrease in slope. This implies 
a slightly nonlinear relationship between and the amount of center fiber 
debonding. Using a second order least squares equation to approximate the 
equation of this line we find 

o = 1.171*1(T 6 + (4.783 x 1(T 8 )P - (3.079xl(T 9 )P 2 6.9 
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6.5) High Temperature Moduli with Center Fiber Debonding 

Center fiber debonding at high temperature is now to be considered. The 

model was changed slightly to allow for symmetric debonding and the constituent 

material properties were adjusted to reflect a use temperature of 1500°F. In 

figure 5.5 the degradation of longitudinal and transverse moduli is represented. 

E is seen to decrease by 7.3% when the center fiber is completely debonded. 
in 

This is less of a decrease than observed for room temperature E , but this 

in 

difference is not significant. The least squares linear approximating equation for 
the line depicting the decrease of E^ in figure 5.5 follows 

E = 4.681xl0 7 - (3.096 x10 5 )P 6.10 

The longitudinal modulus of the fiber alone is over 10 times greater than that of 
the matrix. Since 11.1% of the structural effect of the fiber mass has been 
removed due to the total debonding of the center fiber yet only a 7.3% decrease 
in E is witnessed and considering E is ten times E then the conclusion is 
that the change in E^ brought about by center fiber debonding is significant, 
but not as significant as the change seen in a^. 

E and E are shown to be coincident in figure 5.5. Ignoring geometric 

122 133 

effects and considering that E is 11 times greater than E^, it would be 
expected that as the center fiber becomes structurally inactive and the matrix 
material gains a larger relative percentage in the determination of the 
composite’s properties, there would be an increase in the composite’s transverse 
modulus. Due to geometric effects, the opposite occurs. E is shown to have 
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decreased by 2% when the center fiber is debonded. This decrease is not 
significant in magnitude and is mainly due to the reduced effective 
cross-sectional area, and the subsequent reduced transverse stiffness of the 
composite, as resulting from the debonding of the center fiber. The least squares 
equation describing the decrease is given by 

E ss E = 4.063 x10 6 -(7.879x10 3 )P 6.11 

122 133 V 1 

This decrease may be due to the fact that when the fiber debonds there is, in 
effect, a cylindrical void where the fiber use to be. This hole removes some of the 
stiffness from the model resulting in a lower transverse modulus. Note that this 
decrease is almost of the same percentage as the room temperature decrease in 
E and E . 

122 133 

6.6) High Temperature Poisson’s Ratio with Center Fiber Debonding 

In figure 5.6 the Poisson’s ratios of the composite are plotted. Note again, 
v runs coincident with v , and v is coincident with v . v and v both 

112 113 123 132 112 123 

decrease by approximately 1.5% when the center fiber is debonded fully. The 
least squares linear equation for ^ and is given by 

1 / = 1 / « 0.289 -(3.726x10- 4 )P 6.12 

112 113 \ . l 
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and for v and u it is 

123 132 


!/=!/= 0.256 - (3.718*10‘ 4 )P 6.13 

123 132 


When compared to the room temperature results, it is evident that there is a 
small increase in the change of v and v , and the change in v and v is 

112 113 123 132 

about the same. From these results it may be concluded that temperature has 
little effect on the rate of degradation of these material properties even though 
the actual value of the properties may be different due to the use temperature. 
v , v ,v and v may all be considered to be relatively insensitive to center 

112’ 113’ 123 132 J J 

fiber debonding. 


6.7) High Temperature Shear Moduli with Center Fiber Debonding 


The effect of high temperature fiber debonding on shear moduli is shown 
in figure 5.7. G and G have degraded by almost 5% with the debonding of 
the center fiber. It is important to note that G is 3.5 times greater than G 
and G . In this case, G is 3.5 times greater than G , yet G decreases as 
the matrix material’s structural role increases. The least squares approximation 
of G and G is given by 

112 113 6 J 


G = G = 2.042x10 6 -(8.532x10 3 )P 6.14 

112 113 v ' 
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The change in G and G at room temperature was 3% and at high 
temperature it was 5% when the center fiber was debonded; These changes 
represent the slope of the lines on their respective graphs. The difference in these 
slopes is not significant. 

G is also represented by a line on figure 5.7. G , at high temperature, 

123 123 

had degraded by almost 3% when the center fiber had totally debonded. G is 

m 

5.5 times greater than G . Geometric effects dominate and the results indicate 
' a decrease in G . Again, G at high temperature has been shown to degrade 
by approximately the same amount as at room temperature. 

G = 1.651x10* - (4.093* 10 3 )P 6.15 

123 

G , G and G can all be said to be relatively insensitive to center 

112 113 123 

fiber debonding. 

6.8) High Temperature Thermal Expansion Coefficients with Center Fiber 
Debonding 

Figure 5.8 contains the plots that represent the changes in thermal 
expansion coefficients at high temperature as the center fiber is debonded. Note 
that the lines representing a and a are coincident, a and a increased 

122 133 122 133 

by less than 1%, an insignificant amount. The least squares approximation for 

a and a is given by 

122 133 6 J 
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a = a = 2.174xl0‘ 5 + (1.761*1(T 8 )P 6.16 

122 133 

In figure 5.8, a shows an increase by 17.3% when the center fiber is 
completely debonded which represents a decrease in effective fiber area of 11.1%. 
This substantial increase is due to the fact that a is negative and a is 

fll m 

positive. The least squares approximation of a is given by 

a = 1.709*1(T 6 + (2.612xlCT 8 )P 6.17 

in 

This change in a is significant and close to the percent change in a at room 
in m 

temperature. 

6.9) Room Temperature Moduli with All Fibers Debonding 

The case when all fibers are debonded at room temperature is now 

considered. The effect of all the fibers debonding on the moduli of the composite 

is plotted in figure 5.9. The degradation of E is now obviously non-linear. 

E , with 100% fiber debonding, has decreased by 78.2% and is lower than E . 
in m 

The cylindrical holes left by the debonded fibers decrease the stiffness of the 
composite, much like a sponge. The composite can actually have a lower 
modulus than that of the matrix alone. Though the fibers are still defined in the 
finite element model the connectivity has been removed, in effect producing a 
large cylindrical voids along the debonded length. Displacements may then cause 
the fiber and matrix materials to overlap. Of course this is not realistic and is a 
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result of the way the modeling was done. This affect, may cause the calculation 
of material properties to be to conservative. The least squares quadratic 
approximation of the equation for is given by 


E =4.132*10 7 — (7.685x 10 5 )P + (4.573xl0 3 )P 2 


6.18 


E and E are coincident and appear to be linear as seen in figure 5.9. 

1 x o *5 

E and E decreased by only 15% with 100% fiber debonding. Note that E is 

122 1-33 IP 

19 times greater than E . Therefore the decrease in E and E is due to the 

fil 122 133 

cylindrical voids in the model. The linear approximation for E and E is 

122 133 

given by 


E = E = 7.216*10 6 -(1.069x10 4 )P 

1 2 2 1 3 3 
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6.10) Room Temperature Poisson’s Ratio With All Fibers Debonding 

At low levels of fiber debonding, up until about 25% debonding, v 
increases. Since u is an enforced displacement and does not change, then v must 
be increasing in this range. 
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After 25% debonding, v decreases with increasing debonding. The initial rise 
may be due to the fact that the Poisson’s ratio of the fiber is smaller than that of 
the matrix and E is also less than E . Therefore, as the fibers are released 

f22 m 

there is reduced composite stiffness in the transverse direction. As more fiber is 
released, geometric effects begin to control. Deflections transverse to the fiber 
direction are seen more in the shrinking diameters of the hollow cylinders, 
created by the removal of the continuity between fiber and matrix, than in 
displacements on the external surface of the model where the average 
displacement was calculated. A cubic least squares approximation was used to 
describe this effect. 

v =u = 0.328 + (1.901*10‘ 3 )P 6.21 

- (6.007*10 _5 )P 2 + ( 3.826* 10' 7 )P 3 

v and v are again coincident as seen in figure 5.10. 100% debonding 

123 132 

brought about a decrease in v and v of 26.7%. The value at 100% debonding 
is less than that of the matrix alone. This may be due to the spongy effect of the 
debonded composite, where a significant percentage of the transverse 
displacements take place in the cylinders left by the debonded fibers, resulting in 
less transverse strain. The least squares approximation for the curves of u and 
v at room temperature and 0% to 100% debonding is given by 

132 

v =v = 0.244 -(6.532*10- 4 )P 6.22 

123 132 
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6.11) Room Temperature Shear Moduli With All Fibers Debonding 


Figure 5.11 represents the results of debonding on the shear moduli at 
room temperature. G decreased by 24% when the fibers were completely 
debonded and, as with v , its value is less than that of the pure matrix. The 

123 

least squares approximation of the degradation of G is given by 

123 


°123 = 2 - 690xl ° 8 - (6.606*10 3 )P 6.23 


G and G are again coincident as seen in figure 5.11. Both these 

1 lo 

properties decreased by 25% when 100% debonding had been achieved. As 
before, at 100% debonding, these properties are lower in value than that of the 
matrix material alone, presumably due to the sponge effect discussed earlier. The 
linear approximation of this degradation is given by the following equation 


G = G = 3.350xl0 6 — (8.430*10 3 )P 6.24 

112 113 V 1 

6.12) Room Temperature Thermal Expansion Coefficients With All Fibers 
Debonding 

Figure 5.12 shows the change in the composite’s thermal expansion 
coefficients when debonding is taken from 0% to 100% at room temperature, a 
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increases drastically as the debonding increases. At 100% debonding, has 

increased by 557%. The value of a at 100% debonding is virtually identical to 

in 

that of the matrix alone. The least squares equation for the increase in is 
given by 

a = 1.570xl0’ 6 + (8.357 x10‘ 8 )P 6.25 

ill v 

The degradation of a and a is also depicted in figure 5.12. Both a 

1 22 133 122 

and a degraded by only 14% when all fibers were completely debonded. This 

133 

is a very insignificant amount. The significance in these results is that the final 
value for a and a is virtually equal to a . This would tend to imply that 

122 133 m 

the sponge effect discussed earlier has little effect on thermal expansion. The 
equation for this degradation is given by 

a = a = 1.140 x10- 5 -(1.658x10' 8 )P 6.26 

122 133 

It is significant to note that a , a , and a all approach a when the 

° 111 122 133 m 

composite is approaching a quasi— monolithic state. Here, geometric effects seem 
to have little effect. 
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6.13) High Temperature Moduli With All Fibers Debonding 

Consider now high temperature (1500°F) conditions while debonding takes 

place. In figure 5.13 the moduli of the nine cell model are plotted as the percent 

disbond is increased. Note that is non— linear and degrades rapidly during 

initial debonding but tends to stabilize as debonding approaches 100%. This 

curve is similar in shape to that of room temperature except it is shifted 

down the ordinate. The total decrease in E is 85.6% and its value at 100% 

in 

debonding is lower than that of the matrix alone. Again, this is probably due to 
the sponge effect. The least squares quadratic approximation of E is given by 

E = 3.237*10 7 - (7.352*10 S )P+ (4.727*10 3 )P 2 6.27 

The effect of high temperature does not seem to significantly change the rate of 
degradation of the longitudinal modulus due to debonding but it does change the 
actual value of the modulus at any given degree of debonding. 

E and E are shown to coincide in figure 5.13. E and E both 

122 133 122 133 

degrade by 20.6% at 100% fiber disbond. Note also that E is 11 times greater 

D 

than E and E . The least squares approximation of the plot of E and E 

f22 f 33 122 133 

is given by the following equation 

E = E = 4.084xl0 6 -(8.471xl0 3 )P 6.28 

122 133 
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It is of interest to note that the difference between E and E at 100% 

in 122 

disbond may be attributed to the geometry of the "holes" left by the debonding 
fibers. It would be expected that at 100% debonding the composite would 
display isotropic behavior except for the effect of the "holes" left by the debonded 
fibers. The presences of these cylindrical "holes" also allow the composite moduli 
to be lower than that of the matrix alone. 

6.14) High Temperature Poisson’s Ratio With All Fibers Debonding 

Figure 5.14 shows the effect of fiber debonding on the Poison’s ratios at 
high temperature for all fibers debonding. In this figure v and i/ are 

coincident as well as u and v . As with the room temperature results, u 

123 132 112 

increases initially and then begins to decrease at about 25% debonding. The 
total decrease in v and i/ , from 0% debonding to 100% debonding, is 9.1%. 

i/=i/= 0.330 + (2.633*10" 3 )P - (7.972xl0- 5 )P 2 6.29 

112 113 , ' , 

+ (5.041x10' 7 )P 3 

v and v decrease linearly with fiber debonding as seen in figure 5.14. 

123 132 

These material properties decreased by 29.6% when 100% debonding was 
reached. The equation approximating v and v is given by 

123 132 

v =u = 0.255 -(7.522*10' 4 )P 6.30 

123 132 
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6.15) High Temperature Shear Moduli With All Fibers Debonding 


As in figure 5.11, G and G are coincident in figure 5.15. This figure 
shows a linear decrease in G and G with fiber debonding. G and G 

112 113 112 113 

decreased by 35.7% when the composite was at 100% debonding even though G 

m 

is 3.5 times greater than G . This is the same trend seen at room temperature. 
The least squares approximation of the line representing G and G is given by 


G =G = 2.036*10 6 - (7.264x10 3 )P 6.31 

112 113 ■ ' 


The degradation of G , as seen in figure 5.15, is slightly non-linear both 

123 

at high temperature and at room temperature (figure 5.11). G shows a slight 

increase in slope after about 25% debonding. This increase is so slight; that a 

linear approximation of the line was used, as in the room temperature case. The 

linear approximation of G obtained by the least squares method is given by 

123 

G^ = 1.652x10® - (5.690 x10 3 )P 6.32 

The total decrease in G is 34.2% at 100% disbond with G being 5.5 times 
greater than G . Here again, as in many of the cases presented, there is a 

f 23 

higher value of the matrix’s property yet a decrease in the composite’s property 
is witnessed with increasing levels of debonding,. This trend must again be 
attributed to the effect of the cylindrical "holes'’ left by the debonding fiber. 
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6.16) High Temperature Thermal Expansion Coefficients With All Fibers 
Debonding 

Figure 5.16 depicts the change in the composite’s thermal expansion 

coefficients through the range of 0% to 100% debonding at 1500°F. As the 

percent disbond increased, a significant increase in was evident. At 100% 

debonding, a had increased by 799% and was virtually the same as a and 
111 122 

a at 100% debonding. This implies that the cylindrical "holes" do not affect 

13 3 

the isotropic thermal expansion behavior of the completely debonded composite. 
The linear approximation for this increase in is given by the equation 

a = 2.359xl0* 6 + (1.748xlO’ 7 )P 6.33 

in v ' 

and a are also shown in figure 5.16. These values degrade with the 
increasing debonding of fibers. Though a and a decreased by 10.7% in the 

122 133 

range from 0% to 100% fiber debonding, this is not a significant decrease in these 
values. The decrease is described by the equation 

a = a =2.174 x10' 5 ~(2.370x10‘ 8 )P 6.34 

122 133 . \ / 
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CHAPTER 7 
MICROSTRESSES 


Profiles of microstresses through the composite model were plotted to 
allow for the investigation of the stress states in the fiber and the matrix 
materials under different loadings. These microstresses were normalized with 
respect to the equivalent applied load. The loading conditions used are the same 
ones described earlier to find the composite’s material properties. Some of these 
stress plots are described here and additional results are given in the appendix. 
Plots of microstresses were obtained only for loading conditions which produced 
microstresses that were significant with respect to the failure stress of the 
constituent in which they appeared. 

7.1) Loading: <r e 

XX 

Consider figure 7.1, each line represents the microstresses on each of 
the nine faces of grid points in the model (figure 3.3). This particular plot 
represents the normalized constituent microstresses in the 11 direction as a result 
of an enforced displacement in the XX direction when no debonding has been 
modeled. The stresses shown here are plotted along a line running from the 
center of the model, point "O", to point "A" on the edge of the model. Since the 
center fiber is not loaded, a low stress state can be seen in face 9, the unloaded 
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Figure 7.1 - a ■ Normalized Microstresses, 0-A Direction, 0.0% Debonding, a** Loading 



face. There is also an increase in the stress level on the inside edge of the 
surrounding fibers at face 9. This is the result of the additional load they must 
carry before it is transferred to the center fiber. The lines depicting the stress 
states in faces 8 and 9 are coincident in these figures. Since there are eight 
surrounding fibers that carry the additional load, the effect on each is not very 
great. Further inside the composite the stresses are transferred to the fiber 
through the matrix and an equal load sharing between all fibers is achieved. The 
figure shows that the center fiber has reached a nearly full stress state at about 
half the length of the model. 

Figures 7.2 through 7.9 represent the constituent microstresses at 
varying levels of debonding of the center fiber. As debonding progresses, the load 
can not be transferred sufficiently through the depleted fiber— matrix interfacial 
area. The result is a center fiber that is less structurally active. Stress levels in 
the surrounding fibers increase as a result of the debonding, but since there is a 
community of fibers surrounding the debonding fiber, this increase is minimal. 

A slight stress gradient is formed along the length of the surrounding 
fibers as debonding increases. This stress gradient begins to decrease as total 
center fiber debonding is approached. In figure 7.9, total center fiber debonding 
has been reached and the load that was once carried by the center fiber has been 
transferred to the community of surrounding fibers. At this stage the 
longitudinal load gradients have been eliminated and no significant increase in 
the stress levels of the surrounding fibers exists. Figures 7.10 through 7.12 show 
similar results by plotting the normalized constituent microstresses along the line 
segment 0— C. 
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Figure 7.3 — a Normalized Microstresses, O—A Direction, 2.78% Debonding, a e Loading 
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Figure 7.5 — a Normalized Microstresses, 0— A Direction, 5.56% Debonding, a e Loading 
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Figure 7.8 — a Normalized Microstresses, 0-A Direction, 9.72% Debonding, a e Loading 
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Figure 7.10 — a Normalized Microstresses, O-C Direction, 0.0% Debonding, a c Loading 
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Figure 7.12 - cr Normalized Microstresses, 0-C Direction, 11.11% Debonding, a e Loading 
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Figures 7.13 through 7.15 and 7.16 through 7.18 show the shear stresses 
in the 13, a , along directions defined by the line segments O-A and O-C, 

13 

respectively. These shear stresses transfer the load to the debonding fiber and 
are most significant at the interface of the debonding fiber and the matrix. The 
longitudinal shear stresses tend to be highest at the face where the debonding is 
initiated and the fiber begins to be structurally active. Note that when 
completely debonded, the shear stresses drop to virtually zero because there is no 
load transfer to the center fiber. 

Figures 7.19 to 7.21 also display the effect of center fiber debonding on 
longitudinal shear stresses (o^). These results are similar to those discussed 
above. 

7.2) Loading: <r e 
yy 

Figures 7.22 through 7.30 represent the resulting constituent 
microstresses due to an enforced displacement transverse to the fiber direction. 
These graphs show the microstresses along a line defined by the line segment 
0— A as the amount of center fiber debonding increases. In this transverse 
loading condition the matrix material carries most of the load and the fiber is 
strained the most. This is to be expected due to the relatively low transverse 
modulus of the fiber as compared to that of the matrix material. Since the fibers 
are not structurally dominate when loaded transversely, the effect of the 
debonding of the center fiber is not as significant as in the axially loaded case. 
As debonding progresses, stress risers are evident in the matrix material 
bordering on the debonding fiber. Also, no significant increase occurs in the load 
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Figure 7.30 — a 




carried by the neighboring fibers. Note that when the center fiber is completely 

debonded it still shows a small stress level at the first face of grid points (figure 

7.30). This is a result of keeping these points attached to the matrix so as not to 

remove the fiber— matrix connectivity altogether, thus avoiding singularities in 

the finite element analysis. At this level of debonding, the transverse stresses on 

all the other faces of the center fiber are virtually zero. 

Figures 7.31 to 7.33 portray the transverse microstresses in the 0— B 

direction due to a transverse enforced displacement. In these figures the matrix 

material between the fibers contains compressive stresses as a result of the large 
E 

— rr 5 — ratio. This produces a sharp stress gradient through this area. Even 

f 22 

though the stress levels shown here are not large (note scale in figures 7.31 to 

7.33), this large gradient from tension to compression may be significant in the 

failure of the matrix. As center fiber debonding increases, a increase in this 

compressive stress can be seen. At the same time, each debonded face of the 

center fiber approaches a zero transverse stress state. Only small changes in the 

stress levels of the neighboring fibers are witnessed. 

In figures 7.34 to 7.36 the transverse microstresses along the direction 

of the line segment 0— C as resulting from the loading conditions described above 

are plotted at varying levels of center fiber debonding. As in the plots in the 

0-A direction, the stress levels in the fibers are small compared to the stresses in 

the matrix material. Therefore, center fiber debonding has little effect on the 

E 

stress levels of the rest of the composite. Since — g 51 — is large, most of the 

f 22 

transverse straining of the composite takes place in the fiber and most of the load 
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is carried by the matrix. The step— like increases of stress seen in the matrix 
material are a result of the stress field around the more transversely straining 
fibers. 

Figures 7.37 through 7.39 are the graphical representations of the 
constituent microstresses in the 11 direction plotted along a line defined by the 
O-A line segment with varying degrees of center fiber debonding. These 
microstresses are due to a transverse enforced displacement as described earlier. 
Before fiber debonding is initiated, each pair of faces that are an equal distance 
from the center of the model have the same stress state, as demonstrated by the 
coincident lines for faces 1 and 9, 2 and 8, 3 and 7, and 4 and 6 (figure 7.37). 
These faces are referred to as sister faces. Face 5 is in the plane of symmetry at 
the center of the model and, therefore, has no sister face. In the 11 direction, the 
fiber is in compression and the matrix is in tension. This is due to the large 
E E 

- £ m ratio, the small— g* 1 - - ratio, and the difference in the Poison’s ratios of the 

f 22 fll 

fiber and the matrix. As fiber debonding begins the sister faces no longer have 
the same stress profiles. With increasing fiber debonding, stress levels in the 
center fiber begin to drop as well as those in the neighboring matrix material of 
the debonded faces. For faces where the matrix is still bonded to the center fiber, 
nominal stress risers can be seen forming on the edge of the neighboring fiber and 
in the matrix areas nearest the debonded fiber (figure 7.38). When the center 
fiber is totally debonded, these edges of the matrix and the neighboring fibers 
show signs of local stress concentrations due to the absence of the center fiber 
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(figure 7.39). This increase in stress is not distributed across the area of the 
matrix or the area of the fiber but is manifested in the form of local stress risers 
in each. 

Normalized constituent microstresses in the 11 direction are again 
displayed in figures 7.40 through 7.48 but in these figures the stresses are plotted 
along a direction defined by line segment O-B. Again, these stresses are a result 
of a transverse enforced displacement. Initially, all stresses are compressive and 
sister faces, as described above, maintain consistent stress profiles (figure 7.40). 
With the first layer of elements debonded (1.39% fiber debonding), a sharp stress 
riser appears in tension in face 9 and the stress profiles of the sister faces diverge 
(figure 7.41). As the next layer of elements is allowed to debond, the stress riser 
in the center fiber is diminished. As before, stress risers develop at the edges of 
the neighboring fibers and in the matrix. The stress concentrations in the 
neighboring fibers dominate (figure 7.42). When the fiber is completely 
debonded, the stress level in the center fiber approaches zero and the resulting 
stress concentrations are seen as local effects at the edges of the fiber and matrix 
(figure 7.48). 

Figures 7.49 through 7.51 depict the stress {a ^ ) profiles along a line 
defined by O-C due to the load condition discussed above with varying levels of 
center fiber debonding. Again, for the case of no fiber debonding, stress profiles 
of each pair of sister faces closely correlate with each other. With increasing 
center fiber debonding, the stress levels of the surrounding matrix and fiber 
materials are not significantly affected. Some drop in the tensile stress level in 
the neighboring matrix material can be seen along with some increase in the 
compressive stress level of the neighboring fiber, neither of which are very 
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dramatic. When the center fiber is at 100% disbond the stresses in the center 

* 

fiber have diminished but no significant change has taken place in the stress 
profiles of the rest of the composite. 

Figures 7.52 to 7.54 depict the microstresses in the 33 direction (a ) 

3 3 

due to a transverse load, a . They are plotted along the 0— B direction for 

yy 

various degrees of center fiber debonding. The most dramatic level of stress is at 
the interface between fiber and matrix. As the center fiber is debonded, these 
spikes of compressive stress increase as the stress state in the center fiber is 
relaxed. The drop in the stress state of the center fiber is minimal. Virtually no 
increase in the stress level of the neighboring fiber is seen, even when the center 
fiber is totally debonded (figure 7.54). The result of the debonding is a localized 
increase in compressive stress in the matrix material bordering on the debonded 
fiber. 

Figures 7.55, 7.56, and 7.57 depict the transverse shear microstresses, 
a , developed in the composite as the result of a transverse enforced 

23 

displacement. These microstresses are plotted along lines in the 0— A, 0— B, and 
0— C directions, respectively. No significant change in the stress profiles shown 
in these figures takes place with the onset of center fiber debonding. The 
majority of the shear stresses occur in the matrix regions, therefore, the 
debonding of the center fiber does not significantly effect the stress levels in the 
rest of the composite. 
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7.3) Loading: <y e 

yx 

In figures 7.58 through 7.60 the constituent microstresses ( 0 ^) due to 

an enforced displacement in the YX direction are plotted for increasing levels of 

center fiber debonding along a line in the direction defined by the line segment 

0— A. The matrix portion of the composite is stressed to a greater degree than 

G f G 

the fiber. This is expected because of the large — ^r 2 — ratio s 6 

f 12 L f 12 

each layer of elements in the center fiber is released from the matrix, the stress 
level in that fiber layer approaches zero and an increase in the stress level of the 
matrix material at that point is evident. Little change in the stress profile of the 
neighboring fiber can be seen, even at total center fiber debonding (figure 7.60). 
The change in the stress state, due to the debonding of the center fiber, manifests 
itself in the form of a local shear stress riser in the matrix material at the edge of 
the debonded fiber with little evidence of the increased stress level being 
distributed to the community of surrounding fibers. This is evident when figure 
7.58 is compared to figure 7.60. 

Figures 7.61 to 7.63 are shown to reflect the effect of center fiber 

debonding on constituent microstresses (a ) along a line in the 0— B direction. 

12 

These stresses result from an enforced shear displacement in the YX. Center 
fiber debonding produces the most significant changes in the stress profile in the 
matrix at the point of debonding and, as expected, in the center fiber. These 
areas decrease in their level of stress as the center fiber is increasingly debonded. 
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Some stress relaxation can be seen in the fiber that is not being debonded, though 
this change in the fiber’s stress state is quite small (note scale in figures 7.61 
through 7.63). 

By the use of figures 7.64 to 7.66 an attempt has been made to show 
the effect of center fiber debonding on a microstresses. These figures represent 
the microstresses occurring when an enforced shear displacement in the YX has 
been applied to the composite model. They are plotted along a line defined by 
the line segment 0— C. Center fiber debonding causes the stress levels to drop on 
faces of the center fiber that have been debonded. On these same grid faces, the 
stresses in the matrix next to the debonding fiber decrease, but not as 
dramatically as those in the center fiber. Stress profiles of the neighboring fiber 
remain virtually unchanged as the percent disbond increases. Little shear stress 
sharing takes place among fibers in the composite because most of the cr shear 
stress is carried by the matrix. The decrease in the stress level created by the 
debonding of the center fiber is manifested as local stress reductions at the edge 
of the matrix material nearest the debonding fiber. 

7.4) Loading: <7 e 
yz 

Figures 7.67 to 7.69 depict the effect of center fiber debonding on the 

a microstress profiles, plotted in the 0— A direction, resulting from an applied 
23 

shear in the YZ direction. The cr shear stress is mostly carried by the matrix 

23 

material in the P100— Graphite/Copper composite. As debonding is modeled, the 
stress level of the center fiber decreases as each face of grids in the fiber is 
released from the matrix. The stress level in the matrix material local to the 
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Figure 7.69 - a Normalized Microstresses, O-A Direction, 11.11% Debonding, a e Loading 
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debonding interface also drops, but not as significantly as in the center fiber. 
Little change can be seen in the stress- profiles of the grid faces in the neighboring 
fiber. This suggests that, since the matrix is structurally dominate in shear, little 
change in the shear stress state of the other fibers will be evident. 

In figures 7.70 to 7.72 the same load and microstresses are considered, 
but here they are plotted along a line defined by the O-B line segment. The 
trend here is very similar to that stated above except that the stress levels in the 
matrix material are not as great as before. Debonding produces a slight 
reduction in the local matrix stress state, and little change is seen in the stress 
profile of the other fiber. 

Figures 7.73 through 7.75 are a microstress plots in the 0— C 

23 

direction resulting from a YZ shear at increasing levels of center fiber debonding. 
As before, only a slight decrease in the stress level of the interfacial matrix 
material can be seen as the result of the center fiber debonding. Virtually no 
change in the stress level of the other fiber can be seen as center fiber debonding 
is propagated. Again, the dominate shear modulus of the matrix masks the 
effects of the debonding from the other fibers. 

7.5) Loading: Temperature Load 

The effect of center fiber debonding on longitudinal thermal 
microstresses (cr^) along a line in the direction of the 0— A line segment is given 
in figures 7.76 through 7.84. Figure 7.76 shows the stress profile before fiber 
debonding. Note the effect of the large a mismatch between fiber and matrix. 


149 



D / 


Ot 

o 


EFFECT OF 0.0% FIBER LENGTH DEBONDING 
ON CONSTITUENT MICROSTRESSES «7 23 > 
DUE TO A LOAD IN THE YZ-DIRECTION 




2 




0 


-0.5 


F1BBR 


KW n'saa 


T 


FIBER 


ImmI 


0.000 


0.009 0.010 

O-B AXIS (INCHES) 


0.019 


0.020 


FACES 

t > 1 i t / I , 



I 2 I * 5 4 7 • 
LAVERS 


FACE 1 
FACE 2 
FAC E 3 
FACE 4 
FACE 5 
FACE 6 
FACE T 
FACE8 
FACB 9 


Figure 7.70 — a Normalized Microstresses, O—B Direction, 0.0% Debonding, a e Loading 

23 yz 



D / 


Cn 


EFFECT OF 5.56% FIBER LENGTH DEBONDING 
ON CONSTITUENT MICROSTRESSES ((7 u ) 
DUE TO A LOAD IN THE YZ-DIRECTION 



P 

|Q: 

Q 


IM 


% 

o 


0 

O 


N 

>• 




0.5 


-0.5 


0.000 


F1DBR 


EZME9 


FIBER 


0.00B 0.010 

O-B AXIS (INCHES) 


— i — 
0.018 


3ML 


^LOZO 



LAYERS 

FACE 1 
FACE 2 
FACE 3 
FACE4 
FACE 5 
FACE 6 
FACE T 
FACES 
FACE 9 


Figure 7.71 — a ^ Normalized Microstresses, 0— B Direction, 5.56% Debonding, a e Loading 



/ cr 


Cn 

to 


EFFECT OF 11.11% FIBER LENGTH DEBONDING 
ON CONSTITUENT MICROSTRESSES (CT 23 ) 
DUE TO A LOAD IN THE YZ-DIRECTION 


2-1 





LAVERS 

FACE 1 
FACE 2 
FACE 3 
FACE 4 
FACE 5 
FACE 6 


-0.5 4 


F1BBR 


IZttiillM 


FIBER 


r 


0.000 


0.005 0.010 

n-R AXIS (INCHES) 


— I — 
0.016 


lm\ 


~q!o20 


FACET 
FACES 
FAC E 9 


Fieure 7 72 - a Normalized Microstresses, O-B Direction, 11.11% Debonding, a* Loading 
6 23 y z 



/ cr 


2.5 


EFFECT OF 0.0% FIBER LENGTH DEBONDING 
ON CONSTITUENT MICROSTRESSES ((7 23 ) 
DUE TO A LOAD IN THE YZ-DIRECTION 


t 


* i.s 


rt 

<N 


0.5 




V 


uz 


•0.5 
0.000 


FIBER 


0.008 


MATRIX 

r 


T 


FIBER 


T 


I MATRIX 1 


0.010 0.010 0.020 
O-C AXIS (INCHES) 


0.020 


Toso 



FACE 1 
FACE 2 
FACE 3 
F ACB 4 
FACE 5 
FACB 6 
FACET 
FACB8 
FACB 9 


Figure 7.73 — a Normalized Microstresses. O— C Direction, 0.0% Debonding, a e Loading 
° 23 ys 
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Figure 7.74 - a Normalized Microstresses, O-C Direction, 5.56% Debonding, c e Loading 
° 23 yz 
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Figure 7.75 — <x^ Normalized Microstresses, O-C Direction, 11.11% Debonding, a e Loading 
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Figure 7.77 — a ^ Normalized Microstresses, O— A Direction, 1.39% Debonding, Thermal Loading 
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Figure 7.78 - a Normalized Microstresses, O-A Direction, 2.78% Debonding, Thermal Loading 
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Figure 7.79 - Normalized Microstresses, O-A Direction, 4.17% Debonding, Thermal Loading 
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Figure 7.80 — a Normalized Microstresses, 0— A Direction, 5.56% Debonding, Thermal Loading 
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Figure 7.81 - a Normalized Microstresses, O-A Direction, 6.94% Debonding, Thermal Loading 
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Figure 7.82 — Normalized Microstresses, 0— A Direction, 8.33% Debonding, Thermal Loading 
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Figure 7.83 - a ^ Normalized Microstresses, 0-A Direction, 9.72% Debonding, Thermal Loading 



do/ 


EFFECT OF 11.11% FIBER LENGTH DEBONDING 
ON CONSTITUENT MICROSTRESSES (Cr u ) 2 

DUE TO AN APPLIED THERMAL LOAD |q|q|0] c 

n Qi0- 

000 


FACES 

j i » s i r ■ 


- 260 -f — 
0.000 


MATRIX 


0.010 

0-A AXIS (INCHES) 


I 2 1 » 5 I T • 
LAVERS 


FACE 1 
— FACE 3 
FACE 3 
FACE 4 
FACE 5 
FACE 6 
FACET 

) FACES 

0 020 FACE 9 


Figure 7.84 — a .Normalized Microstresses, 0— A Direction, 11.11% Debonding, Thermal Loading 



The matrix is in compression while the fiber is in tension. This produces large 
stress gradients at the fiber-matrix interface. In this thermal loading case, the 
composite is fixed in the center so it displaces symmetrically. Therefore, each 
line describing the stress state of a given face is coincident with the line 
describing the stress state in its sister face. Fiber debonding is not implemented 
symmetrically. Therefore, lines describing the stress state of each sister face 
diverge with debonding of the center fiber (figure 7.77). 

As debonding progresses, as seen in figures 7.78 through 7.83, stresses 
in the center fiber decrease, particularly those in the middle of the fiber. Stresses 
in the neighboring fibers increase, particularly in the area nearest the debonding 
fiber. When the center fiber is totally debonded, all longitudinal stresses in the 
fiber are virtually eliminated and the load it was carrying has been transferred to 
the surrounding fibers (figure 7.84). Consider that a similar event took place in 
the debonding of the axial loaded composite. The surrounding fibers inherited 
the center fiber’s load when it was debonded. In the thermal loaded case the 
resulting increase in stress in the neighboring fibers is seen more as a local stress 
concentration at the edge of the fibers nearest the center fiber. In the case of the 
axially loaded composite the added load is distributed evenly across the entire 
area of the neighboring fibers (figure 7.9). This difference can be attributed to 
the large difference in the thermal expansion coefficients of the fiber and the 
matrix. Figures 7.85 to 7.87 depict similar results but along the direction defined 
by the line segment O-C. In figure 7.87 the effect of fiber debonding on local 
stress increases in the neighboring fiber is not as significant due to the greater 
distance from the debonding fiber. 
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Figure 7.85 — a Normalized Microstresses, 0-C Direction, 0.0% Debonding, Thermal Loading 
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Normalized Microstresses, O-C Direction, 5.56% Debonding, Thermal Loading 
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Figure 7.87 — a Normalized Microstresses, 0— C Direction, 11.11% Debonding, Thermal Loading 
n 



Figures 7.88 to 7.90 and 7.91 to 7.93 depict the effect of fijDer debonding 
on longitudinal shear microstresses (cr^) due to a thermal load and are plotted in 
the directions defined by the line segments 0— A and 0— C, respectively. These 
figures provided a graphic representation of the high shear stresses between fiber 
and matrix due to the large difference in thermal expansion coefficients of the 
fiber and the matrix. The areas of high shear stress are at the interface between 
the fiber and matrix materials which corresponds to the location of the large 
gradient of axial stresses mentioned earlier. This suggests load transfer between 
fiber and matrix by means of shear stresses. With increasing fiber disbond, the 
stresses in the center fiber are relaxed and the center fiber becomes structurally 
inactive. 

Figures 7.94 through 7.96 also show the effect of fiber debonding on 
longitudinal shear microstresses ( a ) due to a thermal load and are plotted in the 
directions defined by the line segment 0-A. These results are similar to those 
described above. 

It is of interest to note that any loading case which caused the 
cylindrical "holes" left by the debonding fibers to decrease in diameter would in 
actuality cause a compressive stress on the disbonded fiber. This compressive 
stress would induce friction between fiber and matrix and increase the global 
stiffness of the composite. In turn, these effects would change the stress profiles 
discussed above. 
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Figure 7.88 - a Normalized Microstresses. O-A Direction, 0.0% Debonding, Thermal Loading 
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EFFECT OF 5.56% FIBER LENGTH DEBONDING 
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Figure 7.89 — a Normalized Microstresses, 0— A Direction, 5.56% Debonding, Thermal Loading 
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Figure 7.91 - Normalized Microstresses, O-C Direction, 0.0% Debonding, Thermal Loading 
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Figure 7.92 — a Normalized Microstresses, 0— C Direction, 5.56% Debonding, Thermal Loading 
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Figure 7.93 — a Normalized Microstresses, 0— C Direction, 11.11% Debonding, Thermal Loading 
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Normalized Microstresses, 0— A Direction, 5.56% Debonding, Thermal Loading 
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Figure 7.96 — a Normalized Microstresses, 0— A Direction, 11.11% Debonding, Thermal Loading 




CHAPTER 8 
CONCLUSIONS 


From the results of the material property degradation analysis it can be 

seen that, in general, a rule of mixtures based on active fiber volume is not a 

viable method in determining the composite’s material properties. Of the 

material properties considered, the longitudinal thermal expansion coefficient 

seems to be the most sensitive to fiber debonding. This sensitivity of a makes 

in 

it a good indicator of the level of debonding in a composite. If is above 
predicted values then the suspicion of the presence of debonding may be justified 
and the degree of debonding may be approximated by one of the plots in this 
report. 

For small amounts of debonding composite material properties did not 
change significantly. High temperature results indicate a shift in the value of the 
material properties while maintaining the same rate of degradation as in the 
room temperature results. 

When all the fibers were debonded, resulting in a higher percent 
disbond, more significant changes in the material’s properties were observed. 
The longitudinal modulus, E^, and v both show nonlinear degradation at 
these higher levels of fiber debonding. Again, high temperature results show a 
shift in the value of the material properties while maintaining the same rate of 
degradation as in the room temperature results. 
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In the plots of the normalized microstresses the effects of debonding the 
center fiber are considered. In general, the changes in the internal stress state of 
the composite are manifested in two ways: 1) by stress risers local to the 

debonding area and 2) by the neighboring fibers sharing the increased level of 
stress. Varied amounts of each of these effects can be seen in the stress plots. 
The amount of each is dependant on the loading conditions and the composite’s 
constituent material properties. An example of the local stress riser effect is seen 
in the a microstresses generated by an applied YX displacement (figures 7.58 - 
7.60). The distribution of the additional stress to the other fibers can be seen in 
figures 7.1 through 7.9. These figures depict the microstresses in the 11 direction 
produced by a longitudinal enforced displacement. Other examples are discussed 
in the previous chapter. 

In cases where the microstresses in the fibers are small, fiber debonding 
had little effect. An example would be the a ■ microstresses resulting from a 
displacement transverse to the fiber direction (figures 7.55 — 7.57). 
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CHAPTER 9 
SUMMARY 


This work has shown that the greatest effect of fiber debonding on 
composite material properties occurs when neighboring fibers are debonded and 
all load must be carried by the matrix material. This method of finding material 
properties of debonded composites is a viable one yet the amount of required 
computer time is very costly. Some of the extra cost in computer time can be 
attributed to the use of superelements. While they are convenient to use in a 
modular finite element mesh such as this, they are also computer intensive and 
can lead to high analysis cost. The ability to do restarts using the saved 
database did alleviate some computer overhead on the runs where only the 
boundary conditions changed and not the mesh itself. 

Some of the stress plots have large jumps in stress due to the presence 
of a ring of elements with large aspect ratios around each fiber. Where this does 
not significantly effect the material property calculations, it may appear 
conspicuous when considering the stress states in the fiber and the matrix. 

To achieve the results described here, over 230 large MSC NASTRAN 
superelement computer runs were made. This alone is stifling for the average 
computer installation. Many supporting FORTRAN programs were written and 
utilized to aid in the data reduction, material property calculations, 
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least-squares approximations, microstress normalization, and plotting. Each of 
these programs consume additional computer time. This inordinate amount of 
computer utilization suggests that, while yielding consistent results, this form of 
investigation into fiber debonding is very time consuming and expensive. 

The modeling of debonding of composite materials using the finite 
element method provides the ability to obtain the stresses internal to the 
composite. These microstresses can not be attained by experimental methods. 
This method provides unique insight into the mechanics of load transfer between 
fiber and matrix. 
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APPENDIX 


The following pages contain additional microstress plots. These 
provide the ability to analyze the internal stresses state at each level of fiber 
debonding. These normalized microstress plots have been include to complete 
the descriptions of the stress states described in the body of this paper. 
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Figure A.2 - a Normalized Microstresses, O-C Direction, 1.39% Debonding, a* Loading 
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Figure A. 3 — a Normalized Microstresses, O-C Direction, 2.78% Debonding, a* Loading 
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Figure A.5 - a Normalized Microstresses, O-C Direction, 5.56% Debonding, a e Loading 
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Figure A.6 - a Normalized Microstresses, O-C Direction, 6.94% Debonding, a e Loading 
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Figure A.7 - a Normalized Microstresses, O-C Direction, 8.33% Debonding, a 
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Figure A.8 — a Normalized Microstresses, 0— C Direction, 9.72% Debonding, Loading 
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Figure A. 9 — a Normalized Microstresses. O— C Direction, 11.11% Debonding, a e Loading 
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Figure A. 10 — a Normalized Microstresses, O— A Direction, 0.0% Debonding, cr* Loading 
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Figure A. 11 - a Normalized Microstresses, O-A Direction, 1.39% Debonding, a e Loading 
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Figure A.12 - a Normalized Microstresses, O-A Direction, 2.78% Debonding, a e Loading 

° 31 xx 



/ a 


CO 

CO 


EFFECT OF 4.17% FIBER LENGTH DEBONDING 
ON CONSTITUENT MICROSTRESSES ((J 31 ) 
DUE TO A LOAD IN THE XX-DIRECTION 




FACE 1_ 
FACE 2 
FACE 3 
FACE 4 
FACE 5 
FACE 6 
FACB 7 

FACES 
0.020 FACB 9 


Figure A. 13 — 0 ^ Normalized Microstresses, O— A Direction, 4.17% Debonding, a e Loading 




/ cr 


EFFECT OF 5.56Z FIBER LENGTH DEBONDING 
ON CONSTITUENT MICROSTRESSES ((7 31 ) 
DUE TO A LOAD IN THE XX-DIRECTION 



0.2 H 



- 0.1 


-0.2- 


FIBER 


eemi 


T 


T 


FIBER 


0.000 


0.008 0.010 

O-A AXIS (INCHES) 


— I — 
0.018 


ML 


n 

0.020 



LAYERS 

FACE 1, 

FACE 2 
FACE 3 
FACB4 
FACE 5 
FACE 6 
FACE 7 
FACE8 
FACE 9 


Figure A. 14 - a 
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Figure A. 15 - Normalized Microstresses, O-A Direction, 6.94% Debonding, Loading 
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Figure A. 16 — a 
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Figure A. 17 - a Normalized Microstresses, O-A Direction, 9.72% Debonding, a e Loading 
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Figure A.18 - a Normalized Microstresses, O-A Direction, 11.11% Debonding, a e Loading 
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Figure A. 19 - Normalized Microstresses, O-C Direction, 0.0% Debonding, a e Loading 
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.20 - a Normalized Microstresses, O-C Direction, 1.39% Debonding, a e Loading 
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Figure A.21 — a 
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Figure A. 22 — a Normalized Microstresses, O— C Direction, 4.17% Debonding, a e Loading 
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Figure A.23 — a Normalized Microstresses, O-C Direction, 5.56% Debonding, Loading 
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Figure A. 24 — a Normalized Microstresses, O— C Direction, 6.94% Debonding, a e Loading 
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Figure A. 25 — a Normalized Microstresses. 0— C Direction, 8.33% Debonding, a e Loading 
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Figure A. 26 — a Normalized Microstresses, 0— C Direction, 9.72% Debonding, o e Loading 
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Figure A.27 — a Normalized Microstresses, 0— C Direction, 11.11% Debonding, a e Loading 
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Figure A.28 - a Normalized Microstresses, O-C Direction, 0.0% Debonding, <r e Loading 
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Figure A.29 — <r Normalized Microstresses, O-C Direction, 1.39% Debonding, Loading 
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Figure A. 31 — a ^ Normalized Microstresses, O— C Direction, 4.17% Debonding, Loading 
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Figure A.32 - a Normalized Microstresses, O-C Direction, 5.56% Debonding, a e Loading 
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Figure A.33 — a Normalized Microstresses, O— C Direction, 6.94% Debonding, a e Loading 
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Figure A. 35 - a Normalized Microstresses, O-C Direction, 9.72% Debonding, a e Loading 
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Figure A.36 — a Normalized Microstresses, 0— C Direction, 11.11% Debonding, cr e Loading 
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Figure A.37 - a Normalized Microstresses, O-B Direction, 0.0% Debonding, a e 
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Figure A.38 - a Normalized Microstresses, O-B Direction, 1.39% Debonding, a e 
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Figure A.56 — a Normalized Microstresses, O-A Direction, 1.39% Debonding, a e Loading 
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Figure A. 66 - a Normalized Microstresses, O-C Direction, 2.78% Debonding, <r® Loading 
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Figure A. 70 - a Normalized Microstresses, O-C Direction, 8.33% Debonding, a e Loading 
° u yy 



257 


EFFECT OF 9.72% FIBER LENGTH DEBONDING 
ON CONSTITUENT MICROSTRESSES «T U ) 
DUE TO A LOAD IN THE YY-DIRECTION 




I 11*54/1 
LAVERS 


FACB 1 
FACE 2 
FACES 
FACE 4 , 
FACB S 
FACB 6 
FACE 7 
FACBJB 
FACB 9 
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Figure A. 75 — a Normalized Microstresses. O-B Direction, 2.78% Debonding, <r e Loading 
33 yy 
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Figure A.76 - a Normalized Microstresses, O-B Direction, 4.17% Debonding, a Loading 
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Figure A. 77 — a Normalized Microstresses, 0— B Direction, 5.56% Debonding, a e Loading 
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Figure A. 78 - a Normalized Microstresses, 0— B Direction, 6.94% Debonding, a e Loading 
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Figure A.80 — a Normalized Microstresses. O-B Direction, 9.72% Debonding, a e Loading 
33 yy 
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Figure A. 81 - a Normalized Microstresses, O-B Direction, 11.11% Debonding, <r e Loading 
° 33 yy 
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Figure A. 82 — a Normalized Microstresses, O-A Direction, 0.0% Debonding, a® Loading 
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Figure A. 83 — a Normalized Microstresses, O-A Direction, 1.39% Debonding, <r e Loading 

12 yx 
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Figure A.84 - a Normalized Microstresses, O-A Direction, 2.78% Debonding, a e Loading 
5 12 y* 
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Figure A.85 — a Normalized Microstresses. O-A Direction, 4.17% Debonding, <r e Loading 
° 12 yx 
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Figure A.86 - a Normalized Microstresses, 0-A Direction, 5.56% Debonding, a e Loading 
° 12 
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Figure A.87 — a Normalized Microstresses, O-A Direction, 6.94% Debonding, a e Loading 

5 12 yx 
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Fieure A 88 - a Normalized Microstresses, O-A Direction, 8.33% Debonding, a e Loading 
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Figure A.89 - a Normalized Microstresses, O-A Direction, 9.72% Debonding, a e Loading 
° 12 y* 
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Figure A.90 — a Normalized Microstresses, 0— A Direction, 11.11% Debonding, a e Loading 

12 yx 
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Figure A. 91 — a ^ Normalized Microstresses, O-B Direction, 0.0% Debonding, a e Loading 
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Figure A.92 - a Normalized Microstresses, O-B Direction, 1.39% Debonding, a e Loading 
° 12 y* 
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Figure A. 93 — a Normalized Microstresses, 0— B Direction, 2.78% Debonding, a e Loading 
12 yx 
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Figure A. 94 - a Normalized Microstresses, O-B Direction, 4.17% Debonding, a e Loading 
12 yx 
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Figure A. 95 — a Normalized Microstresses, O-B Direction, 5.56% Debonding, a e Loading 
12 yx 
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Figure A.96 - a Normalized Microstresses, O-B Direction, 6.94% Debonding, a Loading 
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Figure A. 97 — a Normalized Microstresses, 0— B Direction, 8.33% Debonding, a* Loading 

° 12 yx 
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.98 - a Normalized Microstresses, O-B Direction, 9.72% Debonding, a e Loading 
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EFFECT OF 11.11% FIBER LENGTH DEBONDING 
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Figure A.99 - a ^ Normalized Microstresses, O-B Direction, 11.11% Debonding, a e Loading 
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Figure A.100 - a Normalized Microstresses, O-C Direction, 0.0% Debonding, a** Loading 
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Figure A. 101 — a 
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Figure A. 102 — a Normalized Microstresses, O-C Direction, 2.78% Debonding, cr e Loading 
12 y* 
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Figure A. 103 — a Normalized Microstresses, O— C Direction, 4.17% Debonding, a e Loading 
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Figure A. 104 - a Normalized Microstresses,. O-C Direction, 5.56% Debonding, a e Loading 
° 12 yx 
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Figure A. 105 — a Normalized Microstresses, O— C Direction, 6.94% Debonding, <r e Loading 
12 y x 
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Figure A. 106 - a Normalized Microstresses, O-C Direction, 8.33% Debonding, a e Loading 
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Figure A. 108 — a Normalized Microstresses, 0— C Direction, 11.11% Debonding, cr e Loading 

° 12 y* 
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.110 — a Normalized Microstresses, O-A Direction, 1.39% Debonding, a e Loading 
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Figure A. 112 — a Normalized Microstresses, O-A Direction, 4.17% Debonding, cr e Loading 
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Figure A. 113 — a Normalized Microstresses, 0— A Direction, 5.56% Debonding, a e Loading 
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Figure A.I14 - a Normalized Microstresses. O-A Direction, 6.94% Debonding, a e Loading 
° 23 y* 
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Figure A. 115 — a Normalized Microstresses, O-A Direction, 8.33% Debonding, a e Loading 
23 y* 
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Figure A.116 — a Normalized Microstresses, O—A Direction, 9.72% Debonding, a e Loading 
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Figure A. 117 — a Normalized Microstresses, O— A Direction, 11.11% Debonding, cr e Loading 
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Figure A. 118 - a Normalized Microstresses, O-B Direction, 0.0% Debonding, a e Loading 
° 23 ya 
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Figure A. 119 - a Normalized Microstresses, O-B Direction, 1.39% Debonding, a e Loading 
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Figure A. 120 — a Normalized Microstresses. O-B Direction, 2.78% Debonding, a e Loading 
** 23 y* 
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Figure A. 121 — a Normalized Microstresses, 0— B Direction, 4.17% Debonding, cr e Loading 
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Figure A. 122 — a Normalized Microstresses, 0— B Direction, 5.56% Debonding, <r e Loading 
b 23 y* 
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Figure A. 123 - a Normalized Microstresses, O-B Direction, 6.94% Debonding, a e Loading 
° 23 y* 
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Figure A. 124 - a Normalized Microstresses. O-B Direction, 8.33% Debonding, a e Loading 
b 23 y* 
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Figure A. 125 — a Normalized Microstresses, O—B Direction, 9.72% Debonding, <r e Loading 
° 23 y* 
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Figure A. 126 — a Normalized Microstresses. 0— B Direction, 11.11% Debonding, cr e Loading 
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Figure A. 127 — a Normalized Microstresses, O— C Direction, 0.0% Debonding, a e Loading 
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Figure A. 128 — a Normalized Microstresses, O-C Direction, 1.39% Debonding, a e Loading 
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Figure A. 129 — a Normalized Microstresses, 0— C Direction, 2.78% Debonding, <r e Loading 
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Figure A.130 - a Normalized Microstresses, O-C Direction, 4.17% Debonding, a e Loading 
° 23 y* 
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Figure A. 131 — a Normalized Microstresses, 0-C Direction, 5.56% Debonding, a e Loading 
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Figure A. 132 - a Normalized Microstresses, O-C Direction, 6.94% Debonding, <r e Loading 
° 23 y* 
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Figure A. 133 — a Normalized Microstresses, 0-C Direction, 8.33% Debonding, <r e Loading 
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Figure A. 134 - a Normalized Microstresses, O-C Direction, 9.72% Debonding, a e Loading 
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Figure A. 135 - a Normalized Microstresses, O-C Direction, 11.11% Debonding, a e Loading 
b 23 yz 
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Figure A. 136 - <7 Normalized Microstresses, O-C Direction, 0.0% Debonding, Thermal Loading 
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Figure A. 137 — Normalized Microstresses, 0— C Direction, 1.39% Debonding, Thermal Loading 



CO 

to 


EFFECT OF 2.7 87 . FIBER LENGTH DEBONDING 
ON CONSTITUENT MICROSTRESSES (CT U ) 
DUE TO AN APPLIED THERMAL LOAD 



QJ 

O 

O 



o 


IQ 


,c 

B Y 


FACES 

2 ) 1 5 4 7 I 


12)1 5 4 7 1 
LAVERS 


FACE 1 
FACE 2 
FACE 3 
FACE 4 
FACE 5 
FACE 6 
FACE IT 
FACE 8 


Figure A. 138 — a Normalized Microstresses, 0— C Direction, 2.78% Debonding, Thermal Loading 
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Figure A. 139 - Normalized Microstresses, O-C Direction, 4.17% Debonding, Thermal Loading 
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Figure A. 140 — a Normalized Microstresses, 0— C Direction, 5.56% Debonding, Thermal Loading 
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Normalized Microstresses, O-C Direction, 6.94% Debonding, Thermal Loading 
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Figure A. 142 — a Normalized Microstresses, O-C Direction, 8.33% Debonding, Thermal Loading 
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Figure A. 143 — Normalized Microstresses, 0— C Direction, 9.72% Debonding, Thermal Loading 
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Figure A. 144 — Normalized Microstresses, O-C Direction, 11.11% Debonding, Thermal Loading 
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Figure A. 145 — a Normalized Microstresses, O— A Direction, 0.0% Debonding, Thermal Loading 
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Figure A. 146 — a Normalized Microstresses, 0— A Direction, 1.39% Debonding, Thermal Loading 
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Figure A. 147 — a 
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Normalized Microstresses, O-A Direction, 2.78% Debonding, Thermal Loading 
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Figure A. 148 — .<r 


Normalized Microstresses, O— A Direction, 4.17% Debonding, Thermal Loading 
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Figure A. 149 — a Normalized Microstresses, O— A Direction, 5.56% Debonding, Thermal Loading 
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Figure A. 150 - a Normalized Microstresses, O-A Direction, 6.94% Debonding, Thermal Loading 
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Figure A. 151 — a 
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Normalized Microstresses, 0— A Direction, 8.33% Debonding, Thermal Loading 



EFFECT OF 9.72% FIBER LENGTH DEBONDING 
ON CONSTITUENT MICROSTRESSES CO* 3l ) 
DUE TO AN APPLIED THERMAL LOAD 


U. 

O 


n 

b 



M 
O 

o 


me 

DO 


FACES 

? 3 1 S * f • 


0.000 


0.005 0.010 0.016 

O-A AXIS (INCHES) 


FACE 1 
FACE 2 
FACE 3 
FACE4 
FACE S 
F ACE 6 
FACET 
FACEjB 
0.020 PACE 9 


LAVERS 


Figure A. 152 — a Normalized Microstresses, O—A Direction, 9.72% Debonding, Thermal Loading 
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Figure A. 153 — a Normalized Microstresses, O-A Direction, 11.11% Debonding, Thermal Loading 
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Figure A. 154 — a Normalized Microstresses, 0— C Direction, 0.0% Debonding, Thermal Loading 
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Figure A. 155 — a 


31 


Normalized Microstresses, 0— C Direction, 1.39% Debonding, Thermal Loading 
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Figure A. 156 — a Normalized Microstresses, 0— C Direction, 2.78% Debonding, Thermal Loading 
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Figure A. 157 — a 
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Figure A. 158 — a Normalized Microstresses, 0— C Direction, 5.56% Debonding, Thermal Loading 
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Figure A. 159 
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Normalized Microstresses, O-C Direction, 6.94% Debonding, Thermal Loading 



do t 


EFFECT OF 8.33% FIBER LENGTH DEBONDING 
ON CONSTITUENT MICROSTRESSES ((J 3 i) 
DUE TO AN APPLIED THERMAL LOAD 


□cp 

000 

ooo 


h r — 1 

t * I i 


r-V 


rK"" 


-50 +- 
0.000 


MATRIX 


0.010 0.016 0.020 
O-C AXIS (INCHES) 


FACES 

i I 1 5 t T • 


I 2 9 1 « t t • 
LAVERS 


FACB 1 
FACE 2 
FACES 
FACB 4 
FACE 5 
FACB 6 
FACB 7 

i FACES 

o.oao pACB 9 


Figure A. 160 — a Normalized Microstresses, 0— C Direction, 8.33% Debonding, Thermal Loading 
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Figure A. 161 — a 
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Normalized Microstresses, O— C Direction, 9.72% Debonding, Thermal Loading 
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Figure A. 162 - a Normalized Microstresses, O-C Direction, 11.11% Debonding, Thermal Loading 
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Figure A. 163 - a Normalized Microstresses, 0— C Direction, 0.0% Debonding, Thermal Loading 
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Figure A. 164 — a Normalized Microstresses, O-C Direction, 1.39% Debonding, Thermal Loading 
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Figure A. 165 - a Normalized Microstresses, O-C Direction, 2.78% Debonding, Thermal Loading 
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Normalized Microstresses, O-C Direction, 4.17% Debonding, Thermal Loading 
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Figure A. 167 a ^ Normalized Microstresses, 0— C Direction, 5.56% Debonding, Thermal Loading 
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Figure A.168 - a ^ Normalized Microstresses, O-C Direction, 6.94% Debonding, Thermal Loading 
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Figure A.169 — a Normalized Microstresses, 0— C Direction, 8.33% Debonding, Thermal Loading 
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Figure A. 170 — a Normalized Microstresses, 0— C Direction, 9.72% Debonding, Thermal Loading 
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Figure A.171 



Normalized Microstresses, 0— C Direction, 11.11% Debonding, Thermal Loading 
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